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I.  INTRODUCTION 


1.1  The  Concept  of  Locking 

Solids  exhibit  wide  and  significant  variations  in  physical  properties. 
Hius  it  is  extremely  difficult  to  characterize  the  mechanical  behaviour 
of  solids  by  a  single  equation  or  a  single  set  of  equations  which  can  be 
used  in  practical  analysis  to  determine  stresses,  strains,  displacements 
and  other  related  quantities.  Therefore  an  engineer  is  compelled  to 
devise  'ideal  materials'  such  as  the  perfectly  elastic  solid,  the  per¬ 
fectly  plastic  body  etc.  The  idealization  depends  on  the  particular 
problem  to  be  solved  and  the  accuracy  desired  in  the  final  solution. 

The  locking  material  considered  in  this  report  is  one  of  the  'ideal 
materials ' .  The  phenomenon  of  locking  which  characterizes  this  ideal 
material  can  be  visualized  by  considering  the  following  model. 


(a) 

(*) 
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pig.  1.1. 


A  spring  AB  of  spring  constant  ks  has  a  flexible  wire  attached 
firmly  to  its  ends.  Hie  cross-sectional  area  A  and  the  modulus  of 
elasticity  E  of  the  wire  are  such  that  the  quantity  EA/L  is  very 
large  compared  with  the  spring  constant  k8  of  the  spring.  If  a  tensile 
force  T  is  applied  to  the  ends  of  the  spring,  the  spring  elongates  by 
a  certain  amount  AL.  As  T  is  increased  AL  also  increases.  For  a 
certain  value  of  AL  =  AL  the  wire  becomes  taut  as  indicated  in  Pig. 
1.1b.'  After  the  elongation  has  attained  the  value  AL,  any  attempt  to 
elongate  the  spring  involves  the  elongation  of  both  the  spring  and  the 
wire  together.  Thus  a  considerably  larger  force  is  necessary  to  produce 
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a  given  amount  of  elongation  than  is  necessary  to  produce  the  same  amount 
of  elongation  when  the  displacement  is  less  than  AL.  The  variation  of 
T  with  AL  is  then  as  shown  in  Fig.  1.2 


Fig.  1.2. 


The  system  is  said  to  be  locked  when  AL  =  AL.  If  the  wire  were  infinitely 
rigid,  the  T  -  AL  curve  would  be  the  curve  AED  instead  of  ABC.  This 
type  of  locking  is  called  ideal  locking. 

A  granular  material  consisting  of  grains  of  various  sizes  constitutes 
another  example  of  a  locking  material.  A  simple  model  can  be  constructed 
if  the  material  is  made  of  grains  of  two  sizes  as  indicated  in  Fig.  1.3. 


Fig.  1.3. 


The  larger  grains  are  assumed  to  be  vary  hard  to  deform  compared  with 
the  smaller  grains.  If  the  model  is  compressed  by  a  uniform  pressure 
p  the  volume  decreases.  The  volumetric  strain  is  denoted  by  3e.  For 
a  certain  value  of  3e  =  3e^  the  model  assumes  the  form  as  indicated 
in  Fig.  1.4. 


Fig.  1.4. 
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It  can  be  seen  from  the  figure  that  the  volume  can  be  decreased  only 
very  slightly  beyond  this  state.  However,  shear  strain  can  be  produced 
by  the  application  of  a  shear  force.  (Fig.  1.5) 


Fig.  1.5. 


Hie  stress-strain  diagram  is  then  as  indicated  in  Fig.  1.6. 


Fig.  1.6. 


Shear 

Stress 


The  material  is  said  to  be  locked  when  e  <=  e^- 


1.2  Application  of  the  Idealization 

In  many  cases,  granular  soil  which  behaves  very  similar  to  the 
second  example  considered  in  section  1.1  can  be  idealized  as  a  locking 
material.  Another  example  of  locking  material  would  be  rubber.  It  is 
a  common  experience  that  it  takes  very  little  effort  to  stretch  a  rubber 
band  upto  a  certain  amount.  With  increasing  deformation  it  becomes 
harder  to  stretch  the  rubber.  Hie  material  resists  greater  load  with 
little  deformation.  This  behaviour  fits  very  well  with  the  concept  of 
locking  as  described  earlier.  Stress  problems  in  rubberlike  materials 
can  be  treated  as  problems  in  locking  media. 

These  are  two  examples  in  which  the  material  can  be  idealized  as  a 
locking  material.  However  this  is  not  an  exhaustive  list. 
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1.3  Loading  and  Unloading,  Conservative  Property 

The  spring-wire  model  will  again  be  considered  in  this  section  to 
discuss  the  loading-unloading  properties  of  locking  materials.  If  a 
tensile  force  T  is  gradually  applied  to  the  ends  of  the  spring,  the 
spring  elongates.  The  loading  process  can  be  represented  by  the  line 
ABC  in  the  T  -  AL  diagram  (Fig.  1.7).  If  the  tensile  force  T  is 
now  gradually  reduced  to  zero,  the  unloading  curve  in  the  T  -  AL 
diagram  is  CBA  provided  the  spring  and  the  wire  are  still  elastic; 
i.e.,  the  unloading  curve  follows  the  loading  curve  in  the  reverse 


direction,  i.e.,  the  system  is  conservative.  Systems  which  follow  a  curve 
like  ABCA^  can  easily  be  invented,  but  unloading  along  a  curve  like 
CAg  would  be  in  conflict  with  the  basic  laws  of  mechanics,  since  it 
would  involve  the  creation  of  mechanical  energy.  In  this  report  only 
conservative  materials  will  be  considered. 

1.4  Locking  Condition,  Locking  Surface 

In  the  case  of  a  spring  model  one  can  easily  state  that  locking 
takes  place  when  the  elongation  AL  has  attained  a  value  AL  or  the 
strain  has  reached  a  value  Al/L.  Similarly,  in  the  case  of  the  granular 
material  locking  takes  place  when  the  bulk  strain  3e  has  become  equal 
to  3e j •  Therefore,  in  general,  the  condition  of  locking  depends  on  the 
state  of  strain  in  the  body.  In  an  isotropic  material  the  locking  con¬ 
dition  should  be  independent  of  the  particular  choice  of  the  coordinate 
system.  It  therefore  depends  on  the  three  invariants  J^,  Jg,  J^,  of 
the  strain  tensor  which  are  defined  in  the  following  way  in  terms  of 
principal  strains  €g, 
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>  (1.1) 


Prom  Eq.  (l.l),  the  locking  condition 

0(J1,J2,J3)  =  0  (1.2) 

can  also  be  expressed  in  terms  of  the  principal  strains 

0(«1» e3)  =  0  (1.3) 

This  equation  represents  a  surface  in  the  principal  strain  space  which 
is  called  the  locking  surface.  The  locking  surface  can  also  be  defined 
in  a  general  six-dimensional  strain  space. 

The  “locking  function"  0  is  defined  to  such  away  that  the  strain  state  of 
the  body  Is  elastic  if  0  <  0.  Then  in  ideal  locking  0  is  always  less  than  or  equal  to  Zero, 
ie.,  0  =  0  in  the  locked  regions  of  the  body  while  0  <  0  in  regions  that  are  still  elastic. 
However  in  a  ncn-ideally  locking  body  0  can  be  greater  than  zero. 

1-5  Volumetric  and  Distortional  Locking 

The  locking  condition  can  depend  on  the  three  invariants  J1»J2 
and  Jj  of  the  strain  tensor.  However,  in  isotropic  bodies  (under  small 
deformations)  the  shear  stress  does  not  produce  any  bulk  strain  and  the 
hydrostatic  stress  does  not  produce  any  shear  strain.  Thus  two  different 
types  of  locking  can  be  defined  in  an  isotropic  body. 

(i)  Volumetric,  locking  which  depends  on  the  invariant  of  the 

strain  dilatation. 

(ii)  Distortional  locking  which  depends  on  the  invariants  Jg,J3 
of  the  strain  deviator. 


J1  "  3  (€1  +  €2  +  e3> 

•*2  ~  Z  ^£1  "  £2^  +  (e2  "  e3^  +  ^€3  ~  £1^2^ 
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Again,  due  to  isotropy,  the  stress  and  strain  deviators  are  related  by 
Hooke's  law  in  a  volumetrically  locked  region,  while  hydrostatic  stress 
and  bulk  strain  are  related  by  Hooke's  law  in  distortionally  locked  regions. 
However,  some  regions  could  be  locked  both  volumetrically  and  distortion- 
ally.  In  the  latter  case  there  is  no  further  change  in  strain  if  both 
the  volumetric  and  distortional  locking  are  ideal.  Hie  possible  types  of 
locking  and  the  stress-strain  relationship  in  locked  regions  will  be 
studied  in  the  next  chapter. 

1.6  Uniqueness 

Locking  materials  are  conservative.  Therefore,  a  unique  solution 
for  an  equilibrium  stress  problem  or  a  dynamic  response  problem  can  be 
expected  under  appnqpriate  boundary  and  initial  conditions.  These  will 
be  studied  in  the  next  chapter. 


II.  GENERAL  THEORY 


II. 1  Strain  Energy  In  Locking  Materials 

Strain  energy  per  unit  volume  In  a  locking  material  is  assumed  to 
be  a  unique  function  of  strains.  This  statement  can  be  easily  Justified 
in  the  case  of  a  non- ideally  locking  body  in  which  stresses  depend 
uniquely  on  the  strains  as  discussed  in  section  1.3.  The  statement  can 
also  be  Justified  in  the  case  of  an  ideally  locking  body  where  the 
stresses  can  increase  or  decrease  by  indefinite  amounts  while  the  strains 
remain  constant  because  no  work  is  done  during  such  changes. 

Further,  if  the  assumption  of  small  deformations  is  made,  the  dila¬ 
tation  3e  in  an  isotropic  locking  body  depends  reversibly  and  uniquely 
on  the  hydrostatic  stress  s.  Also  each  component  of  the  strain  deviator 
depends  reversibly  and  uniquely  on  the  corresponding  component  of 
the  stress  deviator •  The  functional  relationship  is  the  same  for  all 
pairs  of  stress  and  strain  deviators,  i.e.. 


6  “  k  (€11  +  €22  +  €33)  =  e(8) 


8  =  3  (all  +  a22  +  °33)  =  S(e) 


€Ij  =  (eij  '  6ije^  =  f(aij) 


°ij  =  ^°ij  "  BiJS^  "  f  (eij) 


(2.1) 


where  6..  is  the  Kronecker  delta. 

lj 

Certain  simplification  in  the  expression  for  strain  energy  can  be 
achieved  by  using  the  preceding  set  of  equations. 


The  increment  of  strain  energy  per  unit  volume  can  be  written  as 
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du 


-E  (s6u  *  °ij)(ae6ij  *  ds;3>  ‘>3  •  i-2.3 

1.3 


Because 


E  -n  ■ 0  ■  I  asn  ■  °  ] 


du  =  3sde  + 


h  >  a '  de 1 

L  u  ij 


(2.2) 


From  Eqs.  (2.1)  and  (2.2),  one  can  see  that  the  strain  energy  can 
he  spilt  Into  two  parts  u^  and  u' .  The  part  uQ  depends  on  the 
dilatation  e  while  the  part  u'  depends  on  the  strain  deviator. 
Furthermore,  the  strain  energy  is  independent  of  the  choice  of  the 
coordinate  system  in  an  isotropic  body.  Therefore,  u'  must  be  a 
function  of  the  invariants  of  the  strain  deviator,  Jg  and  J^- 

Thus  one  can  see  that 


0 

duQ  =  de  =  3sde  (2.3a) 

*  tjj  "a  *  3T  "j  *  I  3  ■  1'2’3  (2'31>) 

II. 2  Relationship  between  the  Components  of  Stress  and  Strain  Tensors 
in  Terms  of  Strain  Energy  Function 

Principal  stresses  (o^,  Og,  a^)  and  principal  strains  (e^,  Eg,  e^) 
will  be  used  in  the  further  analysis.  Then  the  strain  dilatation  is 

e  -  i  (ex  +  e2  +  e3) 
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and  the  strain  deviator  is 


The  three  components  of  the  strain  deviator  are  not  independent  of  each 
other.  One  of  them  can  be  expressed  in  terms  of  the  other  two: 


+ 


Similar  relationships  for  principal  stresses  yield  the  following  result 
O'  -  -  (o'  +  op  (2A) 

Therefore 

du'  =  o|de|  +  Ogd€^  +  °3de3 

*  (2o^  +  op de|  +  (o|  +  2opd€g  (2.5) 


Further, 

J2  *  5  (€i2+  *2  +  €32)  "  J2(€i'  €2} 

*  €i€263  =  J3(€l’  62} 

Hence 


^2 
^2  "3«J 

^3 

"3 


^2 

d6i+3if  d€2 

dj; 

«i + d*2 


Then 


9 


du' 


(*'  „  8»'  ^ 
l33J  sq  *  33J  sq  j  aei 


1 du1  ^g  5u'  00 3  \  .  ,  /- 

*lSj|  3^  2  (a-6) 

Comparing  Eqs.  (2-5)  and  (2.6)  one  can  see  that 

2_t  +  du1  ^2  Su'  bJ2 
2ai  +  °2  =  SjJ  5^  +  3j£  Sq 


+  2o,  du'  ^2  +  au'  ^3 
’l  +  2o2  =  33J 


The  expressions  for  and  in  terms  of  the  components  of  strain 

deviator  are 


=  el2  +  ei2  + 


e:e 


12 


=-e,2c'  -c'e'2 
1  2  12 


(2.6a) 


therefore 


-i  *  »g  -  <*4  *  ‘g>  (*  -  ‘a  *) 
"i *  2»g  ■  <6i  *  2sg>  (tjj  -  6i  sj) 


(2-7) 


One  can  solve  for  o|,  Og  and  hence  from  Eqs.  (2.7)  and  (2.4). 

Thus 


o'  -  e'  1  +  A  /ei2  Pe'c'  -  2e'^  )  ^r 

1  ‘  1  551  3  (€1  2€1  2  2  2  } 


(2.8a) 


-  10  - 


'  62  S3T  *  5  <-8ti2  '  ati'2  *  S22)  Sjj 


(2.8b) 

(2.8c) 


Also,  from  Eq.  (2.3a) 


s  1  dU° 
3  de 


(2.9) 


Then  the  components  of  the  stress  tensor  can  be  calculated  from  the 
following  relationships : 


>  (2.10) 


Equations  (2.8)-(2.10)  show  how  the  stress-strain  law  depends  on  the 
relationships  expressing  the  strain  energy  in  terms  of  strains,  i.e., 
the  strain  energy  function.  The  stress-strain  law  can  then  be  established 
if  the  strain  energy  functions  u(e)  and  the  J^)  are  known. 

Therefore,  the  investigation  in  the  next  few  sections  will  be  to  deter¬ 
mine  the  strain  energy  functions  in  locking  materials. 


°1  =  °i  + 


°2  =  °2  +  8 


03  =  °3  + 


II. 3  Strain  Energy  Function  in  Volumetric  and  Distortional  Locking 

As  discussed  in  Chapter  I,  a  material  locks  when  a  certain  function 
of  principal  strains  called  locking  function  j^e^e^e^)  reaches  a  cer¬ 
tain  given  value.  One  can  define  the  function  0  such  that  for  0  <  0 
the  material  is  elastic,  and  that  0=0  represents  the  state  of  in¬ 
cipient  locking.  Then  in  an  ideally  locking  material  0  >  0  is  impos¬ 
sible.  However  0  could  be  greater  than  zero  in  a  non- ideally  locking 
body. 
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In  general  the  strain  energy  per  unit  volume  uQ  and  the  strain 
energy  of  distortion  per  unit  volume  u'  in  a  locked  region  can  be 
written  as 


uo  =  1  Ke2  +  U0I 


u'  =  2GJg  +  uj 


(2.11) 


In  these  equations  (9/2)k£  and  2GJ^  would  be  the  strain  energy  of 
dilation  and  distortion  if  the  material  were  deforming  as  an  elastic 
material.  They  are  called  elastic  parts  of  the  strain  energy.  uQj(e) 

and  u'(J,i,  J')  are  the  excess  of  strain  energy  over  the  elastic  values. 

I  3 

These  are  called  locked  parts  of  the  strain  energy.  In  further  analysis, 
uj  will  be  assumed  to  be  a  function  of  only. 

As  mentioned  in  Section  I.5>  stress  and  strain  deviators  are  related 
by  Hooke's  law  in  a  volumetrically  locked  region,  while  hydrostatic 
stress  and  dilation  are  related  by  Hooke's  law  in  a  distortionally 
locked  region.  Then  uQ^  is  zero  in  a  distortionally  locked  region 
while  uj  is  zero  in  a  volumetrically  locked  region. 


Similarly,  the  stress  tensor  in  a  locked  region  can  be  split  into 
elastic  and  locked  parts : 


s  =  3Ke  +  Sj 


Jij 


=  2Ge'  +  c 


U  iji 


(2.1 2) 


3Ke  and  2Ge^j  are  the  elastic  parts  of  the  stress  tensor  which  would 
be  the  stresses  if  the  material  were  deforming  elastically,  s^  and 
°ijl  are  excess  stresses  over  the  elastic  values  in  the  locked 
regions.  These  are  called  the  locked  stress  components.  s ^  is  zero 
in  a  distortionally  locked  region  while  o'  is  zero  in  a  volumetrically 
locked  region. 

Now,  let  us  consider  a  body  in  a  locked  state.  Hie  case  of  non-ideal 
volumetric  locking  will  be  studied.  Ideal  volumetric  locking  is  a  limit¬ 
ing  case  of  non-ideal  locking. 
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In  the  volumetrically  locked  state  0  >  0.  Stresses  and  strain 
energy  are  given  by  Eqs.  (2.1l)  and  (2.12)  with  cr^ ^  =  uj  =  0.  If  the 
strain  state  of  the  body  in  the  locked  region  changes  such  that  0 
becomes  equal  to  zero,  the  locked  region  again  becomes  elastic  wherever 
0  has  become  equal  to  zero.  Ihen  the  expressions  for  the  stresses  and 
strain  energy  should  reduce  to  the  elastic  values  as  0  tends  to  zero, 
i«6»  ) 


lim  un  (e)  =  0 

0-O+  01 

lim  s  (e)  =  0 

0  —  ©*•  1 

0+  is  used  to  indicate  that  the  limits  are  taken  from  the  locked  state. 

Uqj  is  a  continuous  function  of  one  variable,  e.  Further,  u^ 
must  tend  to  zero  as  the  locking  function  0  tends  to  zero.  Then  e 
must  be  a  function  of  0.  Also,  as  explained  in  section  1.5  the  volu¬ 
metric  locking  function  0  depends  only  on  the  invariant  ^  =  e.  Thus 

e  =  e(0) 

0  =  0(e) 

In  these  equations  it  has  been  implicitly  assumed  that  e  and  0  are 
interchangeable,  i.e., 

af/0  (2.14) 


(2.13) 


If  e  is  replaced  by  0  in  Eq.  (2.13)  one  obtains 


lim 

0  -0+ 


uQ|(0) 


=  0 


If  it  is  assumed  that 
0  >  0,  one  can  write 


0  ■  0(e)  >  0 


“0/ 


+  . . . 


allows  a  power  series  representation  for 


(2.15) 
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where 


Ot^>  a2,  ...  are  constants. 


The  result  can  be  summarized  in  the  following  way.  Volumetric 
locking  takes  place  when  0  =  0(e)  >0.  In  a  volumetrically  locked 
region  distortion  continues  to  follow  Hooke's  law  while  the  dilation 
must  be  computed  from  the  formulae  that  apply  in  the  locked  regions : 


u'  =  2GJ£ 

U0f  =ai0  +a2^  + 

0  =  0(e)  >  0  ,  M  /  0 


(2.16) 


Similarly,  by  assuming  that  u'  is  a  function  of  J,'  only,  and 

Z  d. 

that  the  distortional  locking  function  0  is  a  function  of  only, 

one  can  derive  the  following  equations  applicable  in  a  distortionally 
locked  region: 


Ke 


2 


u'  =  2GJ£  +  uj 
Ui  =  ^  +  + 


0  =  0( Jg)  >  0  , 


lim  ul  =  0 
0-0+  1 


lim 

0  —0+ 


0 


) 


(2.17) 
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In  concluding  one  can  summarize  the  achievements  in  this  section 
as  follows.  It  has  been  possible  to  express  the  strain  energy  as  a 
function  of  e,  an<*  the  locking  function  0,  in  volumetric  as  well 
as  distortional  locking.  (See  Eqs.  2.1 6,  2.17)  Sfcus  if  the  locking 
function  for  the  material  is  determined  by  experimental  or  other  methods, 
the  strain  energy  function  can  be  determined.  Once  the  strain  energy 
function  is  known  the  stress-strain  relationship  in  the  locked  region 
can  be  established  from  the  formulae  derived  in  section  II. 2. 

In  further  analysis  it  will  be  assumed  that  the  locking  fur  ion 
is  known  and  the  stress-strain  relationship  in  distortional  and  volumetric 
locking  will  be  discussed. 

II. Stress-Strain  Relationship  in  Distortional  Locking 

From  the  discussion  in  section  II. 3,  one  can  write  the  expression 
for  the  strain  energy  function  in  a  distortionally  locking  material  in 
the  following  way.  For  0(Jg)  <  0 


(2.18) 


u'  =  2GJg  +  Px0  +  pgT  +  ... 


u0  "  f  Ke 


u'  =  2GJg 


For  0(j’  )  >  0 


u0  =  §  Ke2 


Then  from  Eqs.  (2.8)-(2.10)  the  stress-strain  relationship  can  be 
obtained.  For  0  <  0  the  strain  energy  function  and  Eqs.  (2.8)-(2.10) 
give  us  the  well-known  Hooke's  law.  However  for  0  >  0  the  following 
stress-Btrain  relationship  can  be  obtained  using  Eqs.  (2.8)-(2.10)  and 
(2.18). 
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s  =  3Ke 

Comparing  these  equations  with  Eqs.  (2.12) 


However,  from  Eqs.  (2.17) 

{°'u' " 0 

Therefore 


=  0 


because  (see  Eqs.  2.17) 


0 


and  ej,  can  not  all  be  zero.  Then 

”1,  •  V  *1  * 


(2.19a) 


1 6 


2 


(2- 19b) 


u2 1 


& 

dJ„  <- 


°j,  -  6s  *  °^2> 


(2.19c) 


Hie  stress-strain  relationship  in  a  distortionally  locked  region  can 
then  be  expressed  by  the  following  relationships 


-  e[  +  s  *  3Ke  +  2Ge}_  +  2pg0  gr  +  O(02) 


a2  =  oj,  +  s  =  3Ke  +  2Ge£  +  2p,,0  ^  +  °(02) 


/  (2-194) 


c,  =  a' 


+  e  =  3Ke  +  2Ge;  +  2$20  ^  +  O(02) 


The  constants  =  0,  S2>  ...  depend  on  the  particular  material. 

Thus  it  has  been  possible  to  express  the  stress-strain  relationship 
in  a  distortionally  locked  region  in  the  form  of  Eq.  (2.19d)  which  con¬ 
tain  infinite  series  on  the  right-hand  sides.  However  these  equations 
can  be  simplified  by  introducing  the  concept  of  stress  increments.  This 
will  be  discussed  in  the  next  section. 


II. 5*  Concept  of  Stress  Increment 

Let  us  consider  a  point  of  the  body  where  0(jg)  has  Just  attained 
the  value  zero.  Then  the  point  is  in  a  state  of  incipient  locking.  The 
stresses  and  strains  at  the  point  in  the  state  of  incipient  locking  are 
denoted  with  a  suffix  I  >  e.g.,  c^,  etc.  Suppose  the  value  of  0 
at  the  point  increases  by  an  arbitrary  small  positive  quantity  0f,  the 
stress  and  strain  at  the  point  increase  from  the  value  at  the  state  of 
incipient  locking.  These  stresses  and  strains  are  denoted  with  a  suffix  f. 
By  using  the  Eq.  (2.19d),  one  can  then  write  the  following  equations 


°1I  =  3Kel  + 


2G€il 


V  ■  +  ao'if  *  %  ‘if  ♦  <«# 
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i.e., 


dif  "  °u  =  3K(ef  ■  ej)  +  2G(€lf  "  €{i) 

+  2Mf  IjJ  €if  +  °(^)  (2.20) 

The  increments  -  c^,  €j^  -  ,  ef  -  e^  are  now  denoted  by  the 

symbols  Acr^,  Ae^,  Ae.  Further  =  0^  -  0j  is  denoted  by  the  symbol 

Then  Eq.  (2.20)  can  be  written  as 

Ao1  =  3KAe  +  2GAe^  +  20^0  €i  + 

Similarly 

Ao2  =  3KAe  +  2GA£^  +  2 0^6  +  0[^2J 

Aa3  =  3KAe  +  2GAe^  +  2^60  ej  +  O[(A0)2] 

The  quantities  A o^,  Ao2  and  A ct3  are  the  stresses  at  the  point 

in  excess  of  the  stresses  at  the  point  at  the  state  of  incipient  locking. 
They  can  be  split  into  two  parts  in  the  following  way: 

Ao.  =  Aa.  +  A  a,  . 

1  le  1/ 


A02  -  io2e  ♦  fio2( 


“s  '  “l.  *  ^3! 


where 

Ac.  *  3KAs  +  2GAe  ’ 
le  1 

Aa2e  =  3KAe  +  2GAe£ 
Ao^  =  3KAe  +  2GAe3 
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are  the  stress  increments  calculated  from  Hooke's  law  and 
Aoli  =  2  %  e[  +  0(^2f)2 

^2,  -  ♦  OUdf 

Ao3(  ‘  ^  ^  6i  *  °<^>2 

are  the  excess  of  the  locked  stresses  over  the  elastic  values.  If  the 
limiting  values  are  considered  as  t#  tends  to  zero  in  these  equations, 
one  obtains  the  following  result. 

d°i/ - (iy  ei 

^  ^  (iy  *2 

d°3 1  "  (fjj)  *3 

Similarly 

do,  =  3Kde  +  2Gde,' 
le  1 

do-  =  3Kde  +  2Gde' 

2e  2 

do^  =  3Kde  +  2Gde^ 

lhe  quantities  do^j  do^,  da3/  are  call“d  locked  stress  increments, 
do^,  dOge,  do^  are  called  elastic  stress  increments.  The  total  stress 
increments  are  then  given  by  the  following  equations : 

(2.23) 


doL  =  dole  ♦  dc^ 


do2  =  do2e  +  do2/ 


do. 


d03e  +  d03i 


(2.21) 


)  (2-22) 
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lhus, using  the  concept  of  stress  increments, it  has  been  possible 
to  obtain  the  simple  relationships  (2.2l)-(2.23)  between  the  stress  in¬ 
crements,  strain  and  strain  increments.  These  increments  have  been 
calculated  at  the  locking  surface ,  i . e . ,  the  surface  0  =  0  in  the 
principal  strain  space.  These  relationships  are  sufficient  to  discuss 
some  interesting  geometrical  properites  pertaining  to  the  locking  sur¬ 
face  and  to  the  problems  of  ideal  dietortional  locking.  These  will  be 
investigated  in  the  next  few  sections. 

II. 6.  Direction  of  the  Locked  Stress  Increment  Relative  to  the  Locking 
Surface _ 

As  discussed  in  section  II. 4,  distortional  locking  takes  place 

when 

0(J£)  =  0 

Now,  J'  is  a  function  of  €',  e'  and  these  in  turn  depend  on  6  ,e_,e  . 
d  -L  d  12  3 

We  may  therefore  write 

0(jp  =  0[J^(€1,€2,€3)] 

and  may  form  partial  derivitives  like 

|g_ 3 


The  definition  of  Jg  as  given  in  section  II. 2,  may  be  rewritten  in 
the  form 

J2^  We3'  =  15  ^2€1  *  €2  ‘  63^  +  ^2€2  ‘  £1  "  €3^ 

+  <2e3  -  €1  -  €2)2] 

from  which 

dJA  , 

35J  =  3  (2ei  ■  e2  ”  €3}  =  €i 


and  hence 


-  20  - 


Then  these  equations  and  Eqs.  (2.2l)  yield  for  the  stress  increments 


or 

dOj  =  202d0  grad  0  (2.24) 

The  vector  grad  0  is  normal  to  the  surface  0=0.  Therefore  it  can 
be  concluded  that  da^  is  normal  to  the  distortional  locking  surface. 
This  is  an  interesting  result  which  is  useful  in  establishing  the 
uniqueness  of  a  solution  in  equilibrium  stress  problems  or  dynamic  res¬ 
ponse  problems  in  an  ideal  locking  medium. 

II. 7*  Magnitude  of  Stress  Increment  and  the  Locking  Surface 

The  locking  function  0f  corresponding  to  the  state  of  strain 
6f  =  ^€lf'  e2f'  €3p  0an  exPresse<*  in  "terms  of  the  locking  function 
corresponding  to  the  state  of  strain  Tj.  =  €g j >  e3I)  by  the 

following  power  series  expansion 
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If  the  state  I  corresponds  to  the  state  0^  =  0  (i.e.,  incipient 

locking),  and  (clf  -  «2f  -  €2I>  e3f  -  are  replaced  by 

Ae  =  (Ae1,  ASg,  Ae3)  ,  Eq.  (2.25)  can  be  written  as 

0f  *  0f  -  0T  =  (grad  0)^  =  ~  •  &  +  0(|ae|2) 


If  one  considers  the  limit  as  Ae  tends  to  zero,  one  obtains 

d0  =  (grad  0)  •  da 

«  =  €I 

From  Eqs .  (2.24)  the  locked  stress  increment  is  then 

da^  =  202[(grad  0)^  _  •*  •  d"e]grad  0  (2.26) 


This  equation  can  also  be  written  aB 


dOj  =  2P2 1 grad  0|^  =  ^  (n-  de)n 


(2.26a) 


where  n  is  the  unit  vector  In  the  principal  strain  space  normal  to  the 
locking  surface  at  "e  =  Then  n  •  d!e  is  the  component  of  increment 
of  strain  normal  to  the  locking  surface.  Therefore  it  can  be  concluded 
that  the  magnitude  of  dc^  is  proportional  to  the  increment  of  strain 
normal  to  the  locking  surface. 

3he  formula  (2.26)  will  now  be  used  in  the  next  section  to  discuss 
ideal  distortional  locking. 
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II. 8.  Case  of  Ideal  Distortions!  Locking 

From  Eq.  (2.26)  Pg  can  be  expressed  as 


0 


2 


Ida,! 

2(grad  0  •  d?)|grad  0 


i.e.,  Pg  is  proportional  to  the  ratio  of  the  magnitude  of  the  increment 
of  the  locked  part  of  the  stress  to  the  increment  of  the  strain  normal  to 
the  locking  surface.  Further,  the  preceding  equation  is  applicable  to 
both  ideal  and  non-ideal  distortional  locking  materials  when  0=0. 

In  case  of  the  ideal  locking  material  stresses  da^  can  increase 
while  the  end  point  of  the  strain  vector  remains  on  the  locking  surface, 
i.e., 

grad  0  .  d€  =  0 

while  |daj  is  not  zero.  Then  Pg  tends  to  infinity.  However  the 
quantity  d\  =  2Pg[grad  0  •  de]  can  tend  to  a  finite  limit.  The  locked 
stress  increment  can  then  be  written  as 


dOj  =  dX  grad  0 

i.e.. 


(2.27) 


(2.27a) 


Suppose  0(Jg)  is  assumed  to  be 


0(Jp  = 


6 


2 

l 


(2.28) 


The  increment  of  the  locked  part  of  the  stress  can  be  calculated 
from  Eq.  (2.27a): 
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do, , 

=  d\€' 

11 

1 

to2< 

= 

a03< 

-  &*'s 

Then  the  complete  Incremental  stress -strain  relationship  corresponding 
to  the  locking  function  given  by  Eq.  (2.28)  is 


do1  =  3Kde  +  2Gde|  +  dXe| 
da2  =  3Kde  +  2Gde^  +  dXe^ 


do3  =  3Kde  +  2Gde^  +  cLXe^ 


(2.29) 


where  d\  is  determined  (as  explained  in  section  II. 11 )  in  the  particular 
problem  being  solved. 

Another  simple  locking  condition  is 


0  =  0(jp  =  0(e^,€p  =  l€i  -  €j!  -  c  =  0 

ijJ  =  1,2/3 


(2.30) 


From  physical  considerations  it  is  obvious  that  the  quantity  |e^  -  c j | 
to  be  used  in  Eq.  (2.30)  corresponds  to  one  of  maximum  difference  be¬ 
tween  two  principal  strain  deviators.  IhiB  can  also  be  interpreted  as 
the  condition  of  maximum  shearing  strain  for  locking.  The  locked  stress 
increments  in  this  case  are 


doii  =  d^^ 

d02/=d^S^  K  -  €j  1 

d03i  =  d*S^ 

The  stress-strain  incremental  relationship  similar  to  Eq.  (2.29)  can 
then  be  obtained. 
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II. 9*  Case  of  Volumetric  Lockli 

In  this  section  the  stress-strain  relationship  in  volumetries lly 
locking  bodies  will  be  studied. 

From  the  results  of  the  section  II. 3  volumetric  locking  takes 
place  when 

0  =  0(e)  =  0 

and  in  the  volumetrically  locked  region 
u'  =  2GJ^ 


25 


»•  * §  M ♦  • 


The  stress-strain  relationship  in  a  volumetries lly  locked  region  is 
then 


°i-2os;  +  3K8*§v1^*- 

"2  ■  S!xi  +  * !  v 3£  *  ••• 


(2.31) 


By  following  a  procedure  similar  to  that  of  section  II. 5,  stress  incre¬ 
ments  can  be  defined.  Thus 


da..  =  do,  +  da, , 
1  le  1/ 


dc2  .  d<,2e  *  d«2( 


d°3  *  d°3e  +  l83/ 


where 


and 


do,  =  2Gdc '  +  3Kde 
le  1 


da.  =  2Gdc '  +  3Kde 
2e  2 


do_  =  2Gdc '  +  3Kde 
3e  3 


da 


1/ 


d02 /  =  d°3i  =  I  V*  lH 


(2.32) 


(2.33) 


(2.33a) 


As  in  the  case  of  distortional  locking,  a  volumetric  locking  sur¬ 
face  in  principal  strain  space  can  be  defined  in  the  following  way. 


0  =  0(e)  =  0(€x,  e2,  *3)  =  0 
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Ey  following  a  procedure  similar  to  that  of  section  II. 6,  the  normality  of 
the  locked  stress  increment  to  the  volumetric  locking  surface  can  be 
established.  Also  the  magnitude  of  the  locked  stress  increment  is 
proportional  to  the  component  of  strain  normal  to  the  volumetric  locking 
surface.  This  result  is  similar  to  that  obtained  in  section  II. 7-  The 
equation  of  locked  stress  increment  is 

d0j  =  20*2 [(grad  0)-de]  grad  0  (2.34) 

where 


11*10.  Particular  Cases  of  Volumetric  Lockir 


A  simple  volumetric  locking  condition  is 
0  =  -(€l  +  e2  +  «3)  - 


(2-35) 


where  is  a  positive  constant.  Here  locking  is  assumed  to  be  in 

compression  only.  Then  from  Eq.  (2.33a) 

=  6a2(de,  de,  de) 

Therefore , 

da1  =  2Gde^  +  (3K  +  6a2)de 
dOg  =  2Gd€g  +  (3K  +  6a2)de 
da3  =  2Gde^  +  (3K  +  6a2)de 


The  hydrostatic  part  of  the  stress  tensor  is  then 
ds  *  (3K  +  6a2)de 

For  a  material  with  the  stress-strain  diagram  of  the  type  of  Fig.  2.1, 
this  equation  can  be  integrated.  Thus 

s  *  (3K  +  60(g)e  +  Sq 

where 
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II. 11.  Problem  of  Equilibrium-Uniqueness 


So  far  the  Investigations  in  the  chapter  have  been  directed 
towards  studying  stress -strain  relationships  in  locking  materials.  In 
this  section  the  formulation  of  an  equilibrium  stress  problem  and  the 
question  of  uniqueness  under  prescribed  boundary  conditions  will  be 
investigated.  Further  discussion  in  the  thesis  will  be  restricted  to 
ideal  locking  materials  only. 


Let  us  consider  a  body  of  locking  material  enclosed  by  a  surface 
S.  Let  the  surface  traction  P^.  be  prescribed  over  a  part  Sp  of  the  surface  S 
and  the  displacements  v^  be  prescribed  over  the  remaining  part 
S  -  Sp  =  Sy.  Depending  on  the  strain  pattern  the  body  will  have  locked 
regions  and  elastic  regions  Vg.  In  particular  cases  either 
or  Ve  could  be  zero.  The  stress  distribution  a  ,  the  strain  pattern 
€p,  and  the  displacements  u^.  in  the  body  corresponding  to  the  given 
boundary  values  should  satisfy  the  equilibrium  conditions,  appropriate 
stress-strain  relationship  in  the  locked  and  elastic  regions  and  kine¬ 
matics.  These  equations  when  referred  to  a  set  of  cartesian  coordinates 
are: 


(i)  Equilibrium  Conditions 


5t 


33T  +  "Sf 


St 


zx 


+  X  =  0 


dr  do  St 

-JSL  +  _Z  +  +  Y  =  0 

dy  dy  dz 


dT 


xz 


dT 


do 


-jE  +  *  +  Z  =  0 

dy  dz 


(ii)  Stress -Strain  Relationship 

(a)  If  ft  <  0,  the  elastic  stress-strain  relationship  should 
be  satisfied.  They  are  is  Lamp's  constant) 


o  =  3A  e  +  2Ge 
x  e  x 

o  *  3A  e  +  2G€ 
ye  y 
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o  ■  3X  e  +  2Ge 
z  e  z 


t  =  2Ge 
xy  xy 


t  =  2Ge 
yz  yz 


t  =  2Ge 
zx  zx 


(b)  If  0=0,  and  the  material  is  an  ideal  distortionally 
locking  material,  the  stress-strain  relationships  are 


do  =  3X  de  +  2Gde  +  d^  . 
x  e  x  de 


do  =  3X  de  +  2Gde  +  dX  , 
ye  y  ae 


do  =  3X  de  +  2Gde  +  dX 


dr  =  2Gde  +  dX 
xy  xy  de 


xy 


dx  =  2Gde  +  dX  5— 1 
yZ  yz  de 


yz 


dT  =  2Gde  +  dl\  , 
zx  zx  de 


zx 


and 


=  0 


(c)  If  0=0,  and  the  material  is  an  ideal  volumetrically 
locking  material,  the  stress-strain  relationships  are: 


do  =  3X  de  +  2Gde  +  dX  , 
x  e  x  ae 


do  =  3X  de  +  2Gdc  +  dX 
ye  y 


-  30 


do. 


It 


xy 


4t 


yz 


dT_ 


zx 


3A  de  +  2Gde 
e  ' 

:  2Gde 

xy 

:  2Gde 

yz 

:  2Gde 

zx 


0  =  0 


(iii)  Kinematics 
£x  = 

V  =¥ 

dw 

z  ~  &z 

V  =  i(|?  +  i ;) 

6yz  =  5  (s  +  S) 
ezx  =  \  (M  +  It) 


From  the  sets  (i),  (ii),  and  (iii)  of  these  equations,  one  can  see  that 
the  equations  and  unknowns  are  well  balanced  in  elastic  and  locked 
regions.  In  an  ideally  locked  region  six  stress-strain  relations,  the 
locking  condition,  six  kinematic  relations,  and  three  equilibrium 
equations  are  available  to  determine  six  stress  components,  six  strain 
components,  three  components  of  displacements  and  the  proportionality 
parameter  A  or  A. 

The  preceding  formulation  of  equilibrium  stress  problem  is  in 
cartesian  coordinates.  Similarly, the  problem  can  be  formulated  in  other 
coordinate  systems.  Now  the  question  of  uniqueness  of  the  solution  of  a 
problem  of  equilibrium  will  be  considered. 
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Statement  of  Uniqueness  Theorem:  Let  a  body  made  of  locking 
material  enclosed  by  a  surface  S  be  subjected  to  the  surface  traction 
Pj  prescribed  over  a  part  8p  of  the  surface  S.  Let  the  displacements 
Vp  be  prescribed  over  the  remaining  part  of  the  surface  S  -  Sp  *  Sy. 
Further,  let  the  stress  distribution,  the  strain  pattern,  and  the  dis¬ 
placements  in  the  body  corresponding  to  these  boundary  values  be  o^, 

€_,  and  u_.  Now,  suppose  the  boundary  values  are  given  increments 
dP  and  dv^,  then  a  unique  stress  distribution  +  do  and  a  unique 
strain  pattern  +  de  are  assured  in  the  body. 

Proof:  Uniqueness  of  +  do,  €p  +  d!c  is  assured  if  the  unique- 
-♦  *♦  «♦ 
ness  of  do,  de  is  proved.  Suppose  the  increments  do,  de  are  not 

unique,  at  least  two  sets  of  solutions 


(i) 

doa  , 

cd 

♦5 

(il) 

dob  , 

are  possible.  Both  solutions  satisfy  the  same  set  of  boundary  values. 
Then  the  stress  distribution 


the  strain  pattern 

de  -  de,  *  e* 
a  b 

the  displacements 
duQ  -  du^  =  u* 

the  solutions  of  a  problem  satisfying  the  following  boundary  values. 

(i)  Surface  traction  on  the  boundary  Sp  is  zero 

(ii)  Displacements  on  the  boundary  are  zero 

Further,  the  body  forces  are  assumed  to  be  prescribed  throughout  the  body. 
Then  they  are  zero  for  the  problem  defined  by  *  quantities. 


and 


are 


(2.37) 
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The  stresses  o*  are  the  solutions  of  a  problem  of  equilibrium. 
Therefore,  they  should  satisfy  the  equilibrium  conditions.  These  con¬ 
ditions  when  referred  to  a  set  of  cartesian  coordinates  can  be  written 
in  the  following  form: 


y1  0*  +  T-  T*  +  T”  T*  =  0 

ox  x  oy  yx  dz  zx 

(2.38a) 

T*  +  T—  0*  +  T*  =0 

ox  xy  oy  y  oz  zy 

(2.38b) 

y-  T*  +  5—  T*  +  T—  0*  =  0 

ox  xz  oy  yz  oz  z 

(2.38c) 

Multiplying  Eq.  (2.38a)  by  u*,  (2.38b)  by  V *,  Eq.  (2.38c)  by  w*, 
adding  and  integrating  the  resulting  expression  over  the  volume,  one 
obtains 


Volume 


rrr  /  do*  dr*  dr*  \ 

JJJ 


dr* 
zx 


+  V 


+  w* 


(- xF'sf'-s* j 

/St*  dr*  da*\ 


dxdydz  =  0 


(2.39) 


Using  the  divergence  theorem  this  equation  can  be  simplified  to  the 
following  form. 


//  [  u*  ( to*  +  m-r*  +  nr*  ) 

JJ  x  xy  xz 


Surface 


Volume 


+  v*  (/t*  +  mo*  +  nT*  ) 

xy  y  yz 

+  w*  (  It*  +  nr*  +  no*)  ]  ds 
zx  yz  z 


[0*€*+  0*6*  +0*6*  +2t*  e*  +2t*  6*  +2t*  6*  ] dxdydz  »  0 
xx  yy  zz  xyxy  yz  yz  zx  zx 


(2.40) 
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where  I,  m,  and  n  are  the  direction  cosines,  and  €*  ...  c*  are 

X  sx 

the  components  of  the  strain  tensor  corresponding  to  the  stresses  a*  .. . 

t*  • 
zx 

Because  the  boundary  values  for  this  problem  are  zero,  Eq.  (2.40) 
simplifies  to 


JJJ 

Volume 


a*  •  e*  dxdydz  «  0 


(2.41) 


In  general  depending  on  the  strain  state,  the  volume  of  the  body  is 
divided  into  elastic  regions  Vg  and  locked  regions  V^.  Then,  Eq.  (2.4l) 
becomes 


a*  •  €*  dxdydz  + 


a*  •  £*  dxdydz  *  0 


V 

e 


V 


I 


(2.42) 


The  positive  definiteness  of  the  integrands  of  this  equation  will  now 
be  studied. 

For  practical  convenience  the  integrand  over  the  elastic  regions 
Vg  is  denoted  by  Ig  and  the  integrand  over  the  locked  region  by 
1^.  Because  the  stresses  and  strains  are  related  by  Hooke's  law  in  the 
elastic  regions 

I  >  0  (2.43) 


In  this  equation  the  equal  sign  applies  only  if  all  ode  are  identi¬ 
cally  zero.  Using  Eqs.  (2.37),  the  integrand  1^  can  be  written  as 

1^  -  (doa  -  do^)  •  (dcg  -  d«b)  (2.45) 

The  type  of  locking  could  be  ideal  or  non- ideal.  Only  the  case  of 
ideal  locking  will  be  considered. 


de  and 
a 

(Pig.  2.2) 


Ideal  locking:  In  case  of  the  ideal  locking  strain  vectors 
de^  in  the  locked  region  must  be  along  the  locking  surface 


The  stress  increments  in  thelocked  region  do  and  da.  can  be  split 

fl  0 

into  elastic  and  locked  stress  increments.  Thus 


da  =  da  +  da  . 
a  ae  at 


da,  =  da,  +  da,  . 
b  be  b  t 


Then  the  expression  for  1^  is 


I,  *  (da  -  da,  )•  (de  -  de,  )  +  (da  -da  )-(d!e  -  dc  ) 
l  ae  be  a  b  at  b  t'  a  b' 


Then,  I.  can  be  written  as 


Xl  *  J1  +  *2 


where 


I..  =  (do  -  da,  ) *  ( de  -  de,  ) 
1  ae  be'  a  b' 


Z2  =  (d‘af  “  ^bi),(d%  ‘  dS} 


(2.46) 


The  stress  (do  -  da,  )  corresponding  to  (de  -  dc,  )  can  be  found 
from  Hooke's  law.  From  the  results  of  section  II. 5>  (^°a/  “  dobj)  18 
normal  to  the  locking  surface  while  (Fig.  2.2)  die  -  de.  is  tangential 

ft  D 

to  the  locking  surface.  Then 
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(2.U7) 


h  ■  (1“„  *  4V'(dI.  •  aS>  '  0 

From  this  equation  and  Eq.  (2.46),  Eq.  (2.42)  can  be  rewritten 
as 


a*  .  e*  dxdydz  + 


o* 

e 


V 

e 


e*  dxdydz  =  0 


where 


a* 

e 


In  V  ,  a*  .  c*  (which  when  expanded  in  cartesian  coordinates  is 

0*€*  +  0*€*  +  o*c*  +  2t*  €*  +  2t*  6*  +  2-i*  €*  )  is  twice  the  strain 

xx  yy  zz  xyxy  yzyz  xz  xz 

energy  2ug  per  unit  volume  of  the  material^  corresponding  to  the  stress 

distribution  a*.  Similarly  a*  •  e*  represents  twice  the  strain  energy 

e 

2u .  in  V.  because  the  contribution  from  the  locked  stresses  are  zero 
le  l 

(see  Eq.  2.47)  and  a*  =  do  -  do  is  related  to  7*  =  di  -  de, 
e  ae  De  a  d 

by  Hooke's  law.  Then  the  preceding  equation  becomes 


2udxdydz  =  0 


where  u  is  the  strain  energy  per  unit  volume  of  the  body.  Then  the 
integral  can  be  zero  only  if  u  is  identically  zero  in  the  body.  Hence, 
only  rigid  body  displacements  u*,  v*,  w*  are  possible.  This  means  that 

o*  3  0 
€*  =  0 


throughout  the  body.  Then 


*Ref.  2.,  p.  171. 
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-♦  -♦ 
doa  =  do^ 

d€  =  de 
a  b 

That  is,  the  increments  do,  d€,  are  unique. 

Though  the  stresses  dc^  do  not  contribute  towards  any  work,  it 
should  be  noted  that  they  are  not  arbitrary  but  are  determined,  through 
the  locking  parameter  dX,  by  the  sets  of  equations  (equilibrium,  stress- 
strain  relationship,  locking  condition  and  kinematics)  as  discussed 
earlier. 

This  concludes  the  proof  of  the  uniqueness  theorem  as  stated  on 
page  32. 


11.12.  Limitations  on  Displacement  Boundary  Conditions  in  Ideal  Locking 

It  has  been  proved  in  the  last  section  that  if  a  solution  of  an 
equilibrium  boundary  value  problem  exists,  it  is  unique.  In  this  section 
a  particular  case  will  be  illustrated  where  the  solution  does  not  exist. 

Consider  a  body  made  of  an  ideally  locking  material.  Let  only 
displacements  be  prescribed  on  the  surface  of  the  body.  Let  the  distri¬ 
bution  of  these  displacements  have  a  fixed  pattern  on  the  surface  and 
increase  the  magnitude  of  v  gradually  from  zero.  For  a  certain  value 

■4  -♦ 

of  v  =  Vj  the  body  locks  in  some  region.  In  general  the  complete  body 
will  not  be  locked  at  once.  As  v  is  further  increased  the  locked  re- 
gion  increases  until  at  a  certain  value  of  v  =  v£,  the  complete  body  is 
locked.  This  is  called  the  'completely  locked  state'  of  the  body.  In 
certain  cases  the  body  can  change  at  once  from  the  elastic  to  the 
completely  locked  state. 

If  the  body  obeys  the  following  distortional  locking  condition 
H  -  C  =  0 

the  body  becomes  completely  locked  when 


throughout  the  body.  Then,  because  the  locking  is  ideal 
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AJg  =  0 


throughout  the  body,  i.e., 


Ae 


3" 


Now  consider  the  strain  energy  u  of  the  body. 


u  =  uQ(e)  +  u'(jp 


After  the  body  is  complete  locked  u'(Jg)  cannot  be  increased,  because 
the  locking  is  ideal.  Therefore  additional  work  can  be  done  only  by 
increasing  uQ(e),  that  is,  the  body  is  deforming  without  distortion. 

But,  the  distribution  of  displacements  on  the  surface  has  a  fixed  pattern. 
Then  the  prescribed  displacement  should  be  such  that  the  body  is  deformed 
without  distortion.  However,  this  is  a  contradiction,  because  such  dis¬ 
placements  will  not  produce  distortional  locking.  Hence  it  can  be  con- 
eluded  that  no  further  increase  of  v  1b  possible  after  the  body  is 
completely  locked  distortionally. 

Similar  results  can  be  obtained  for  volumetric  locking  with  the 
locking  condition  defined  by  the  following  equation 

0(e)  =  0 


11.13.  Dynamic  Problem  -  Uniqueness 

In  the  present  section  the  formulation  of  a  dynamic  response 
problem  and  the  uniqueness  under  certain  boundary  and  initial  conditions 
will  be  studied. 

The  formulation  of  the  dynamic  response  problem  is  very  similar 
to  the  formulation  of  the  equilibrium  problem,  excepting  that  the 
equilibrium  conditions  are  replaced  by  the  equations  of  motion.  The 
second  problem  under  investigation,  the  uniqueness,  will  be  considered 
in  the  following  way. 
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Statement  of  Uniqueness  Theorem  (dynamic  problem);  Consider  a  body 
of  a  locking  material  enclosed  by  a  surface  S.  Let  the  surface  traction 
?(t)  be  prescribed  over  a  part  Sp  of  the  surface  S  and  the  displace¬ 
ments  v(t)  over  the  remaining  part  S  -  Sp  *  Sv«  Also  let  the  stress 
distribution  aj,  the  strain  pattern  £  ,  the  displacements  u 
and  velocity  in  the  body  be  known  at  time  t  =  tp.  Then  during  the 

interval  t  =  tj  and  t  =  tp  +  dt  unique  increments  d?  of  cfp  d? 
of  ?j,  and  du  of  up  are  assured. 

Proof :  If  the  solution  is  not  unique,  at  least  two  solutions  are 
possible. 

(i)  doa,  d€fl,  dufl 

(ii)  d0b,  de^,  dUh 


Both  solutions  correspond  to  the  same  set  of  boundary  values.  Therefore, 
the  stress  distribution 


da  -  da,  =  a* 
a  b 


the  strain  pattern 


dc  -  dc  =  6* 
a  b 


and  the  displacements 

are  the  solutions  of  a  problem  satisfying  the  following  boundary  values. 


du  -  du*.  =  u* 
a 


(2.U8) 


(i)  surface  traction  on  Sp  is  zero 

(ii)  displacements  on  Sv  are  zero 


Furthermore,  body  forces  for  the  problem  defined  by  Eq.  (2.48)  are  zero 
because  the  body  forceB  are  assumed  to  be  prescribed  in  the  body. 
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o*,  €*,  u*  are  the  solutions  of  a  problem  in  dynamics.  Therefore, 
they  should  satisfy  the  equations  of  motion,  appropriate  stress-strain 
relationships  and  kinematics.  These  equations  are  the  same  as  those 
discussed  in  section  11.12,  excepting  that  the  equations  of  motion 
replace  the  equations  of  equilibrium. 

The  equations  of  motion  for  the  stress  distribution  0*  and  the 
displacements  u*  are 


do* 

sr + 


dr* 

T? 


dr* 

,  zx 
+  "sr  - p 


dr*  do* 


drzy  _  d2v* 
“  "  P  dt2 


(2.49) 


dr*  dr*  do* 

ST2  *  T?  *  3T  -  ■= 


d2w* 


By  using  these  equations  an  integral  of  the  following  type  can  be  written 


rrr  law  t  a2».  K  ^  iriz\ 

JJJ  gr 


At 
Volume 


A  dv*  ,  w. 

l#S?" 


/  ^2  _  dT*  do*  dr*  \ 

p^-T?V-aP) 


that  is. 
At 


.  aw  /a2w  th  T1 
♦sr  (s?- ‘ -sr - - s- ]) 


dxdydz  =  0 


At 


M[#  (?♦£♦£) 


j  .  /dr*  do*  dr*  \  v  -  /dT*  dT*  do*\  “1 


dxdydz 


-  40  - 


(2.50) 


fusing  the  divergence  theorem  and  observing  that  the  boundary  values  are 
zero,  this  equation  simplifies  to 

JJJ  (Ix  +  I2)dxdydz  =  0  (2.51) 

Volume 

where 

[(If) **($**  (If)  ] 

J2  =  ’  3t 

I1  can  be  written  as 

-  §  [(iff  *  *  (if)2]  t=t } 

at  t  =  t^  all  the  Increments  u*,  v*,  w*  and  their  time  derivatives 
are  zero  because  we  know  the  solution  cr^,  7^.,  and  ^  at  t  =  t^. 
Then 

Ix>0  (2.52) 

Now  consider  Ig.  We  write 


Further,  at  t  =  t-j- 


d€  =  d€.  =  0 
a  b 
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(2.53) 


therefore , 


(daa  - 


®b>' 


(d?a- 


d€.  ) 

°  tj+At 


The  expression  on  the  right-hand  side  of  this  equation  has  been  dis¬ 
cussed  in  section  11.12.  It  has  been  shown  that  the  integrand  corresponds 
to  twice  the  strain  energy.  By  following  a  procedure  similar  to  that  of 
section  11.12  it  can  be  shown  that  Eq.  (2.5l)  implies  that  only  rigid 
body  displacements,  independent  of  time,  are  possible  and  hence  the 
uniqueness  of  stresses  and  strains. 
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III.  PROBLEMS  OF  BJUILIERIUM  VOLUMETRIC  LOCKING 


HI'l  ^>^-ane  Stress  Problems  In  Volumetric  Locking 

Let  us  consider  a  thin  plate  in  a  state  of  plane  stress  as  shovn 
in  Fig.  3.1.  The  conditions  of  equilibrium  for  the  stresses  are 


Obese  equations  are  identically  satisfied  if  a  stress  function  X  is 
defined  such  that 


o 

x 


A 

7 


a2x  a2x 

°y  =  £5  ’  V  =  *  33y 


(3.2) 


If  volumetric  locking  occurs, 
are  related  in  the  following  way- 


stress  and  strain  in  the  locked  region 


2G<s  - «,) 

2G(«y  -  .,) 


0 

X 


0 


y 


T 

xy 


e 

x 


+  € 

y 


+  6 

Z 


/ 


(3.3) 
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These  equations  can  be  simplified  into  a  more  useful  form: 


a  *  2G(2e  +  e  -  e.) 
x  x  y  l 

a  =  2G(e  +  2e  -  e.) 
y  x  y  l' 


i  e# 

£  =  ts1  ( 2o  -  a  )  +  «— 

x  6G  x  y  3 

€  =  rs  (-0  +  2o  )  +  rr— 

y  6G  x  y  3 


(3.4) 


(3.5) 


Further,  in  order  to  assure  single  valued  displacements  in  a  simply 

connected  body,  the  strains  must  satisfy  the  compatibility  conditions. 

It  is  not  possible  to  satisfy  all  the  compatibility  conditions  if  one 

assumes  that  the  stresses  a  ,  a  and  r  are  functions  of  x  and 

x'  y  xy 

y  only.  It  becomes  necessary  to  assume  that  the  stresses  depend  on 
z  also.  However,  as  in  the  case  of  elastic  materials  ,  it  can  be  shown 
that  the  dependence  on  z  becomes  negligable  as  the  thickness  of  the 
plate  becomes  small.  Then,  the  only  conqiatibility  condition  to  be 
satisfied  by  the  strains  is 


(3.6) 


Then,  by  using  Eqs.  (3.l)  and  (3.4)  this  equation  can  be  written  in 
terms  of  stresses 


v  (o  +  a  )  =  0 
x  y 

where 


V2 


(3-7) 


Ref. 


3,  p.  241. 
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Then,  this  equation  and  Eq.  (3.2)  yield  the  following  differential 
equation  for  X 

V^X  =  0  (3.8) 

Thus,  if  the  stress  function  X  is  obtained  by  solving  this  equation  with 
appropriate  boundary  conditions,  stresses,  strains  and  hence  displacements 
can  be  calculated. 

The  plane  stress  problem  in  an  elastic  locking  material  can  now  be 
stated  in  the  following  way.  In  general,  a  thin  plate  in  a  state  of  plane 
stress  may  contain  locked  and  elastic  regions.  The  stress  function  and 
hence  the  stresses  can  be  obtained  by  solving  Eq.  (3.8).  The  solution  in 
the  elastic  region  can  be  determined  by  methods  of  solving  plane  stress 
problems  in  elastic  bodies.  Then  the  solution  is  complete  if  the  stresses, 
strains  and  displacements  in  elastic  and  locked  regions  satisfy  the  appro¬ 
priate  boundary  and  interface  conditions. 

III. 2  Pure  Bending  of  Beams  in  Plane  Stress 
(i)  Elastic  Solution 

We  shall  now  study  the  stresses  in  a  beam  of  narrow  rectangular 
cross  section  made  of  ideal  volumetric  locking  material  subjected  to  pure 
bending  moments  M  as  shown  in  the  Fig.  3.2.  When  the  value  of  M  is 
very  small  the  corresponding  strains  are  very  small.  The  beam  then  deforms 
as  an  elastic  beam.  The  stresses  in  such  a  beam  are 


a  =o  =  t  —  T  =t  =0 
y  z  yz  zx  xy 
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Fig.  3.2. 


y 


where  h  is  the  depth  of  the  beam  and  b  is  the  width  of  the  beam. 
From  the  preceding  expression  for  the  stresses  one  can  evaluate  the 
dilatation: 


e 

x 


+  e  + 

y 


e 

z 


y 


Because  the  beam  is  made  of  volumetrically  locking  material,  locking 
takes  place  when 


This  locking  condition  assumea  that  the  material  can  lock  in 
compression  as  well  as  in  tension.  Then,  the  beam  subjected  to  pure 
bending  moments  M  locks  at  y  =  i  h/2  when 


M  =  Mf 


6Eh2b 
1  -  2v  t 


(3-9) 


(ii)  Problem  after  Incipient  Locking 

If  M  is  increased  beyond  the  limit  given  by  Eq.  (3-9), 
locked  regions  will  develop  near  the  upper  and  lower  edges  of  the 
cross-section.  There  is  an  elastic  core  between  these  two  locked  regions. 
The  interfaces  between  the  elastic  core  and  the  locked  regions  are  at 
y  =  1  yf  as  shown  in  Fig.  3.3. 
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Fig.  3.3. 


Because  the  locked  regions  are  situated  symmetrically  with  respect  to 
the  x-axis,  the  stress  distribution  is  as  shown  in  the  figure  and  the 
neutral  axis  is  the  line  y  =  0.  The  stress  function  is  governed  by 
the  differential  equation 

V^X  =  0 


in  the  elastic  as  well  as  in  the  locked  regions.  The  boundary  conditions 
for  the  problem  will  now  be  discussed. 


Because  the  stress  and  displacement,  fields  are  antimetric  with 
respect  to  the  neutral  plane  y  =  0,  it  suffices  to  restrict  the  further 
discussion  to  the  region  y  >  0  only.  Then,  (Stresses,  strains,  and 
displacements  in  the  locked  region  will  be  denoted  with  a  bar  over  the 
letters,  e.g.,  ox)  for  x  =  0  and  x  =  L 


T 


xy 


T 

xy 


=  0 


dy  = 


0 


(a) 


(b) 


y 


dy  =  M 


(c) 
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for  y  =  h/2 


T  =0 

xy 


for  y  =  y. 


u  =  u 


v  =  v 


a  =  a 

y  y 


T  =  T 

xy  xy 


€+€+€  = 
x  y  z 


1  -  2V 


o  =  e . 
x  l 


(a) 


(e) 

(f) 

(g) 

(h) 
U) 

(3.10) 


(lii)  Solution  In  the  Locked  Regions 

The  stress  function  X  In  the  locked  region  Is  assumed  as 


X 


a  3 

=  5^ 


Ge 


/  2 
-  y 


(3.11) 


The  stresses  in  the  locked  region  y  >  0  are  then  given  by  the  following 
expressions. 


Hie  corresponding  strains  can  be  calculated  from  Eq.  (3.5)>  Uiey  are 


e 

x 


c 

y 


7 


xy 


=  0 


(3.13) 
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These  three  equations  can  be  integrated  to  obtain  the  displacements. 
Thus, 


“  =  W  *  Wy  + 


2  2  e, 

V=-£2L-S£-  +  -iy+U>x+C 

12G  5T  2  y 


J 


(3.14) 


(iv)  Solution  in  the  Elastic  Region 

In  the  elastic  region  the  stress  function  is 


X 


d  3 


5y 


where  d  is  a  constant.  From  this  the  stresses  in  the  elastic  region 
can  be  calculated: 


°x  *  ** 


o 


y 


(3.15) 


The  corresponding  strains  and  displacements  are  given  by  the  following 
equations : 


€ 

X 


c 

y 


7 


xy 


=  o 


»  .  SH  .  „  y 

E  1' 


+  C, 


2  2 
vdy  dx 

=  "  2E  ‘  2E 


V 


+  c„ 


(3.16) 


(3.17) 
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(v)  Matching  the  Boundary  Conditions 

The  boundary  conditions  (3.10a,b,d,g,h)  are  automatically 
satisfied  by  the  stresses  as  given  in  subsections  (iii)  and  (iv).  lhe 
conditions  (3.10e,f)  yield  the  following  relationships: 


d  a 
1  =  3G 


-Vi 


Cl  =  - 


wy^  + 


(3.18) 


2 

ayl 

ISg 


y  +  C  +  U)X  = 


2. 

It 

2E 


+  C„ 


V 


Equation  (3.19)  can  be  satisfied  only  if  w  =  w^.  Then 


(3.19) 


C 


2 


-  ayl 
'  C2  =  SG" 


From  this  equation  and  Eq.  (3.10i)  we  have 


(1  -  2v) 


yl  ~  &yl  3G 


(1  -  gv) 


C2  '  \  ~  ¥  6  lVl 


(3.20) 


(3.21) 


Now  the  only  condition  to  be  satisfied  is  (3.10c).  This  will  be  dis¬ 
cussed  in  the  next  sub-section. 

From  Eqs.  (3. l8)-(3.20)  one  can  observe  that  three  of  the  five 
quantities  can  be  arbitrarily  chosen.  These  correspond 

to  two  rigid  body  displacements  and  a  rigid  body  rotation  permissible  in 
a  plane  stress  problem. 
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(vi)  Relationship  between  Bending  Moment  and  Curvature 

From  Eqs.  (3.10),  (3.12)  and  (3.15)  the  stresses  in  the  locked 
and  elastic  regions  can  be  expressed  in  the  following  way. 


a 

x 


2L. 

y* 


3Ee , 


o  = 

X 


2(l  +  v)(l  -  2V ) 


Ec, 

L.  _  _± 

2(1  + 


vT 


(3.22) 


By  substituting  these  expressions  in  (3.10c)  one  obtains  the  bending 
moment 


bEc, 

m  =  6TTT7T 


_1 _ hf  3h2") 

1  -  2v  yf  ‘  T  J 


(3.23) 


Strains  in  the  elastic  and  locked  regions  can  be  obtained  using  Hooke's 
law  and  Eq.  (3.5)-  Ihen  ex  in  the  elastic  as  well  as  the  locked  region 
is  given  by  the  following  equation 


2_ 

y* 


(3.24) 


This  equation  shows  that  the  plane  sections  before  bending  remain  plane 
after  bending.  Also,  from  this  equation,  the  curvature  of  beam  is 
found  to  be 


K  = 


1  -  2v 


y* 


Then  the  expression  for  the  bending  moment  becomes 

3 


M 


bE 


6(1  +  V) 


(1  -  2v)  K 


1  3.3  3  ,  ,  2 

1  ~2  +  ¥  h  *  ’  S  e/h 


(3.25) 


variation  of  M  with  k  is  shown  in  Fig.  3.4. 


Pig.  3.4b. 
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The  value  of  M/k  as  tends  to  zero  can  be  evaluated  from 
Eqs.(3.24)  and  (3.25).  It  Is 


M  bEh3 

*  =  8"(T"+  v7 


EI 

e 


(3.26) 


where  I  is  defined  as  the  effective  moment  of  inertia.  For  a  value 
e 

of  V  =  0.25, 


EI 

e 


Ebh3 

10 


(3.27) 


which  is  not  an  appreciable  change  from  the  elastic  flexural  rigidity 

7 

Ebh  /12.  This  can  also  be  seen  from  Fig.  3.4b. 

(vi)  Application  of  the  Results  to  the  Solution  of  Beam  Problems 

As  an  application  of  the  derived  results,  the  problem  of  a 
beam  subjected  to  a  loading  other  than  pure  bending  will  be  considered. 
In  particular,  the  deflections  of  a  simply  supported  beam  AB  under  the 
action  of  a  concentrated  load  P  at  mid  span  will  be  studied  (Fig.  3-5) 


Fig.  3.5-, 
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The  bending  moment  diagram  is  shown  in  the  figure.  According  to  Eq.  (3.9) 
the  beam  locks  when 


M  '  M/  ~  1  -  2v  I 


Then  the  beam  begins  to  lock  at  x  =  L/2  when 


96Eh2b€ . 


P  =  P.  = 


i  “  (1-  2V)L 


(3.28) 


For  P  >  Pj  there  will  be  a  locked  region  xQ  <  x  <  (L/2)  +  Xq.  At  xQ 
the  bending  moment  is  just  equal  to  M^.  Then 


48Eh2be . 


K0  ~  (1  -2v)P 


(3.29) 


If  one  assumes  that  the  moment -curvature  relationship  derived  for  pure 
bending  can  be  used  in  this  case,  i.e. ,  the  change  of  curvature  due  to 
shearing  forces  is  negligible,  one  has  (see  Eq.  3.23) 


x  =  EIk 


0  <  x  <  x„ 


P 

2  x  =  5T 


bE  f  1.3^3  3  „  ,  2 

TTVj  7  1  h  *  ‘  *  e/h  J 


x0  <  x  <  2 


(3.30) 


For  small  deflections  k  can  be  expressed  as  the  second  derivative  of 
the  deflection  w : 


K 


(3.31) 


This  equation  and  Eq.  (3.30)  can  be  used  to  determine  w.  However  as 
explained  in  sub-section  (vi),  even  in  the  extreme  case  where  the  complete 
cross-section  of  the  beam  is  locked,  the  change  of  flexural  rigidity  from 


the  elastic  value  El  is  very  small.  Then  one  can  get  upper  and  lower 
bounds  for  the  deflection  by  considering  the  following  two  cases 

(i)  lhe  complete  beam  is  elastic 

(ii)  Hie  region  0  <  x  <  Xq  is  elastic  while  in  the  region 

Xq  <  x  <  (L/2)  the  beam  is  locked  over  the  entire  cross- 
section. 

In  case  (i)  the  deflection  is 

”  ■  •  si  (x3  - !  l2*)  (3-32> 

In  case  (ii)  we  have  from  Eq..  (3.26) 


Then 
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In  the  extreme  case  when  xQ  =  0 


w  = 


P 

12EI/ 


3 


(3.3 6) 


In  the  two  extreme  cases  given  by  Eqs.  (3.36)  and  (3.32)  the  maximum 
deflections  are  at  x  =  (l/2).  The  values  are 


w 


max 


PL 

55eI 


max 


PL3 

SSeT^ 


For  V  =  0.25,  l/lj  is  given  by  Eq.  (3.23).  It  is  equal  to  0.8333. 

Wien 

w 

=  0.8333 
w 

max 

Die  maximum  deflection  is  reduced  if  the  beam  is  made  of  locking  material. 
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III. 3  Rotating  Disks 

In  this  section  stresses  in  a  rotating  disk  (of  constant  thickness) 
made  of  ideally  locking  material  will  be  studied.  The  present  study  is 
restricted  to  small  deformation. 


(i)  Elastic  Solution 

When  the  angular  velocity  of  the  disk  is  very  small  the  disk 
deforms  as  an  elastic  body.  The  stresses  in  such  a  solid  disk  are* 


a 

r 


a 

z 


3  +  v  2  .2  2v 
— n-i  pai  (R  -  r  ) 


3  +  v  2„2  1  +  3V  22 

pa)  R  -  — -  pa)  r 


Tr0  =  Trz  =  Tez  =  ° 


j 


(3.37) 


where  p  is  the  mass  per  unit  volume  of  the  material,  2R  is  the  diameter 
of  the  disk  and  r  is  the  radial  coordinate.  From  Eqs.  (3.37)  the  dil¬ 
atation  is  found  to  be 


e 

r 


+  +  e 


1  -  2v 

E 


(o 


°0) 


-  2v 


2 

pu) 

E 


[(3  +  v)R2  -  2(1  +  V)r2] 


(3  38) 


Its  maximum  occurs  at  r  =  0  and  is 


(€r  +  €e  +  ez^  =  pu)2r2 

r  ®  z  max  4 


1  -  2v 
E 


It  reaches  the  locking  limit  when 


2 


0) 


4Ee/ 

(3  +  v)pR2(l  -  2 v) 


(3.39) 


Ref. 


3,  p#  TO# 
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(ii)  Problem  after  Initial  Locking 

If  w  is  increased  beyond  Wq,  the  disk  can  no  longer  be 
analysed  as  an  elastic  disk.  Ihere  will  be  a  locked  region  r  <  k  and 
an  elastic  region  k  <  r  <  R.  The  problem  then  reduces  to  that  of  com¬ 
puting  the  stresses,  strains,  and  displacements  in  the  elastic  and  locked 
regions  such  that  they  satisfy  the  equation  of  motion,  the  appropriate 
stress-strain  relationship  and  kinematics  as  well  as  the  boundary  and 
interface  conditions. 


Fig.  3.6. 


(iii)  Expressions  for  Stresses,  Strains  and  Displacements  in  the 
Locked  Region 

Under  the  assumption  of  plane  stress  and  axisymmetry  the 
quantities  to  be  determined  are  the  stresses  ar>aQ>  strains  €r>€e,€z 
and  the  radial  displacement  u. 

Because  the  locking  is  ideal  the  following  locking  condition 


■r*VV  (,r  *  V  '  -e(  lS-k0) 
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holds  throughout  the  locked  region.  In  the  locked  region,  the  stress 
and  strain  deviators  are  still  related  by  Hooke's  law,  i.e., 


2G<Fe  -  €z)  =  a0 
2G(Tr  -  "  °r 


(3.41) 


Other  equations  to  be  satisfied  by  the  stresses,  strains  and  displace¬ 
ments  are  the  equation  of  motion  and  of  kinematics,  i.e., 


da 

r 

IT 


2 

pu  r 


=  0 


e 

r 


du 

dr 


(3.42) 


(3.43) 


Equations  (3.40)-(3.43)  are  sufficient  to  determine  the  unknown  quantities 

V  V  V  V  ?Z  and  "* 

5y  solving  (3.40)  and  (3.4l)  for  stresses  one  obtains 


or  =  2G(2?r  +  e0  -  €t) 
oQ  =  2G(?r  +  2?@  -  tt) 


(3.44) 


When  this  equation  and  Eq.  (3.43)  are  substituted  into  the  equation  of 
motion  a  differential  equation  for  u  results. 


du 


It  has  the  solution 


2  3, 


V) 


(3.45) 
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Then,  the  stresses  and  strains  are  given  by  the  following  expressions : 


r  A  C2  3pw2r2(l  +  v) 
6r  =  C1  '  “5  -  16E - * 

r 


7  „  .  C2  pw2r2(l  +  v) 

*  0,  +  —7:  ■  rrrr"  — ^ 


se  -  ^1  T  “5 

r 


-  E  C2  7o^2r2(l  +  v)  .1 

°r  =  TTT  1_3C1  -  H? - T5t -  '  €/J 

a  -  fsc  +  ^  .  »A2(l  +  v)  .  £  1 

e  1  +  v  L  1  ^  ^L6E  fj 


(3.46) 


(3.47) 


(iv)  Stresses,  Strains  and  Displacements  in  the  Elastic  Region 
From  the  known  solution*  of  the  problem  one  can  write 


‘r  *  C3  -  -5 
r 


€e  =  c3  +  “5 

r 


HT~ 

_  2  2,, 

2s 

3pw  r  (1 

-  v  ) 

8E 

2  2,, 

2\ 

pu  r  (1  - 

JLi 

E  r  E  %  pu2r2(3  +  v) 

1-V  3"l+V  2  "  8 

r 


E  „  .  E  C4  pw2r2(l  +  3v) 

e =  c3 +  rrr  "2  ’  ^ir - 

r 


(v)  Boundary  Conditions 


(3.48) 


The  following  boundary  and  interface  conditions  must  be 
satisfied  by  the  stresses,  strains  and  displacements 


Ref .  3,  p.  70 
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*1 

H 

O 

u  =  0 

(a) 

r  =  k 

u  =  u 

(b) 

r  =  k 

II 

e> 

(c) 

r  r 

r  =  k 

e  +  e.  +  e  =  e. 

(d) 

r  0  z  1 

r  =  R 

■Q 

II 

o 

(e) 

(3.49) 


Hien,  from  Eqs.  (3.45)-(3.48)  and  the  boundary  conditions  (3.49)  one  can 
evaluate  the  constants  C^,  Cg,  and  to  obtain 


Hie  unknown  radius  k  can  be  obtained  from  the  following  equation 
(see  Eq-  3.49c) 


M1  +  y) 

1  -  2V 


2(5  +  V] 
1-  2V) 


pw2  (1  -  2v)2  R2 


(3.51) 


(vi)  Jump  in  oq 

Hie  variations  of  a  and  oa  with  r  for  k/r  *  0.5  is 
r  0 

shown  in  Fig.  3.7b.  In  Fig.  3.7B  pw^^/Ee^  has  been  plotted  against 
k/r  for  several  values  of  -V •  The  purpose  of  the  figure  is  to  study 
the  effect  of  Poisson's  ratio. 
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Fig.  3.7. 
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One  can  see  that  there  is  a  jump  in  Og  at  the  elastic  locking 
interface,  This  can  be  explained  in  the  following  way.  At  the  inter¬ 
face  we  have  satisfied  the  continuity  of  the  displacement  and  the  con¬ 
tinuity  of  the  radial  stress,  i.e.,  for  r  =  k 

u  =  u  a  =  o  (3*52a,b) 

r  r 


From  Eqs.  (3.43),  (3.44) 


E  /  0  du  u  -  \ 

Jr  “  1  +  v  ^  Sr  r  l  J 


(3.53) 


and  from  Hooke's  law 


E  ( du  ,  vu  ) 

'777  7?*“/ 


(3-51*) 


Now,  from  (3. 52b)  °r  =  °r  at  r  =  Then,  from  Eqs.  (3.54),  (3.55) 
and  (3.52a)  one  obtains  at  r  =  k 


3u  du  _  Jr  /,  v  .  J 
3?  "  3?  =  “  2E  (1  *  v  +  v 


!>-l  (t-0 


(3.55) 


This  shows  that  there  is  jump  in  du/dr  at  r  =  k.  Equations  similar 
to  (3.54)  and  (3.55)  can  be  written  for  o@.  Then 


°e  ‘  °e 


uE  1  -  2v  ,  E  /  V  du  du,  €t  \  c^:^ 

=  T  777+rTT  ^TTv3?’3F  +TTv)  {Z'36) 

From  this  equation  and  (3-55)  one  can  see  that  there  is  a  Jump  in  Og. 

(vii)  Problem  After  the  Disk  is  Completely  locked 

As  the  angular  velocity  w  is  increased  the  disk  becomes 
completely  locked  when  (see  Eq.  3.51) 


pu 


1  -  2V 
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After  the  disk  is  completely  locked,  i.e.  for 


1  -2v 


the  stresses,  strains  and  displacements  in  the  disk  are  given  by  Eqs. 
(3.^5)-(3*V7).  Now,  the  boundary  conditions  are 

r  =  0  u  =  0 

r  =  R  a  =  0 

r 


Then,  the  constants  C  and  Cg  can  be  evaluated.  In  this  case  they 
are 


The  expressions  for  the  stresses  can  then  be  written  as 


Similarly,  the  corresponding  strains  and  displacements  can  be  obtained. 
The  variation  of  o  and  aa  are  shown  in  the  Pig.  3.8. 
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IV.  PROBLEMS  OF  EQUILIBRIUM  -  DISTCRTIONAL  DOCKING 


IV. 1  Stresses  Around  a  Small  Spherical  Cavity  in  a  Body  Subjected  to 

Uniform  External  Pressure 

We  shall  consider  a  body  of  distortlonal  locking  material  subjected 
to  uniform  external  pressure  (see  Fig.  k.l)  with  a  spherical  cavity  of 
radius  a.  The  radius  a  is  assumed  to  be  very  small  compared  with  the 
dimensions  of  the  body.  In  the  absence  of  the  spherical  cavity  the 
stress  distribution  in  the  body  corresponds  to  one  of  uniform  hydrostatic 
stress  throughout  the  body.  The  spherical  cavity  modifies  this  stress 
distribution.  However  the  effect  of  the  cavity  is  not  felt  appreciably 
farther  away  from  the  cavity,  i.e.,  at  distances  from  the  center  of  the 
cavity  which  are  large  multiples  of  a.  The  modification  of  the  uniform 
hydrostatic  stress  distribution  in  the  body  due  to  the  spherical  cavity 
will  now  be  investigated. 

(i)  Elastic  Solution 

Let  us  consider  a  large  sphere  of  radius  b  concentric  with 
the  sphere  of  radius  a.  The  radius  b  is  assumed  to  be  very  large 
compared  with  a.  The  modification  of  the  stress  distribution  due  to 
the  presence  of  the  cavity  of  radius  a  is  very  small  at  the  external  surface 
of  the  large  sphere  of  radius  r  =  b.  Then  it  can  be  assumed  that  at 
r  =  b  »  a  the  stress  distribution  corresponds  to  one  of  uniform  hydro¬ 
static  stress. 

As  a  final  step  in  the  solution,  the  derived  formulae  are  modified 
by  letting  the  radius  b  go  to  ».  These  simplified  formulae  are  use¬ 
ful  in  application  to  bodies  whose  dimensions  are  very  large  compared 
with  the  radius  of  the  spherical  cavity  a. 

In  the  elastic  hollow  sphere  of  inner  radius  a  and  outer  radius 
b,  the  stresses,  strains  and  displacements  should  correspond  to  the 
following  boundary  values 


r  =  a  o=0 

r 

r  =  b  a  =  -p 

r 


(4.1) 
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The  geometry  and  the  boundary  values  are  spherically  symmetric.  The 

stresses  a  ,  oa,  the  strains  e  ,  ea  and  the  displacement  u  should 
r  o  r  w 

then  satisfy  the  equilibrium  equation  (oQ>  e^,  are  the  normal  stress 
and  normal  strain  in  any  direction  in  the  plane  normal  to  the  radius 
vector  r) 


do  o  -  o 

_ -  +  2  — - -  =  0 

dr  r 


Hie  stress-strain  relations 


o  +  2a.  =  3K  (c  +  2e.) 
r  0  r  0 


o  -  o  =  2G(e  -  e.) 
r  0  r  0 


(a) 

(b) 


(4.2) 


(4.3) 


and  the  kinematic  relations. 


€  = 


du 

dr 


u 

€0  =  r 


(4.4) 


By  substituting  (4.3)  and  (4.4)  in  the  equilibrium  equation  (4.2)  one 
obtains  the  following  differential  equation  for  u: 


d  u  2  du  2u  » 

— 5  +  7  d?  •  “5  =  0 

dr  r 

Integration  of  Eq.  (4.5)  yields 
C2 

u  =  V  +  -5 
r 


(4.5) 


(4.6) 


Hie  corresponding  stresses  and  strains  can  then  be  calculated  from 
Eqs.  (4.4)  and  (4.5).  Hius, 


2C2 

C - == 

1  ^ 
r 


€e  =  ci  +  "§ 

r 


(4.7a) 
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a  *  3KC,  -  4g 

r  1 


1  "3 

r 


(4.7b) 


°e  =  3KCi  +  20  “5 

r 


From  the  boundary  conditions  (4.l)  and  the  stresses  (4.7b),  one  obtains 


3V3 

C2  =  "  &  “3  3 

b  -  a 


,  =  _IL_b _ 

'1  3K  .3  3 

b  -  a 


(4.8) 


If  one  considers  the  limit  as  b  •*  «,  one  obtains 


C2  =  '  ’  C1  "  ‘  3K 


(4-9) 


'ftien,  one  can  write  the  following  formulae  for  u,  e^,  eQ,  o r  and  oQ 
in  the  body  which  deforms  like  a  completely  elastic  body: 


u  =  .  E-  r  _  EL.  i_ 

3K  4cT  2 


e  =  .  E_  +  1 

r  3K  2G  ^ 
r 


Ik  -  SS“  ^3 


(*•10) 


3 

— -  2L 
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However,  the  body  is  made  of  distortionally  locking  material.  Therefore, 
locking  takes  place  when  the  characteristic  locking  function 

0  =  0(J£)  =  0 


In  this  case  the  locking  condition  is  assumed  to  be  of  the  following  form 

.2 


<‘r  -  *  <‘8  -  +  (e0  ■  •  2£<  ■  0 


(4.11) 


Because  €g  =  e^,  this  equation  simplifies  to 


€r  "  €0  =  1  €l 


(4.12) 


From  Eqs.  (4.10)  and  (4.12)  one  can  see  that  the  body  starts  locking  at 


r  =  a 


when 


4Gc, 


P  = 


(4.13) 


In  the  next  sub-section  the  stresses  in  the  body  after  incipient  locking 
will  be  investigated. 


(ii)  Problem  after  Incipient  Locking 

If  the  pressure  is  increased  from  the  value  given  by  Eq.  (4.13) 
the  locked  region  in  the  body  increases.  Because  of  spherical  symmetry, 
one  can  assume  that  the  region  a  <  r  <  k  (Fig.  4.l)  is  locked  while 
the  region  r  >  k  is  still  elastic,  k  is  the  radius  defining  the 
elastic-locking  interface  when  the  external  pressure  p  >  4Gc  /3.  Then, 

Jv 

the  stresses,  strains  and  displacements  in  the  locked  region  should 
satisfy  Eqs.  (4.14)  through  (4.17),  where  the  barred  letters  denote 
the  quantities  in  the  locked  region. 

(a)  Equilibrium  Equation 


(4.14) 
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(b)  D1 la tat Iona 1  Elasticity 


or  +  20^  =  3K(Ir  +  2€q) 


(U.15) 


(c)  Ideal  Locking  Condition 


€r  "  €e  =  €i 


(4.16) 


(d)  Kinematics 


V! 


7  .$! 

r  dr 


(4.17) 


Die  expressions  for  stresses,  strains  and  displacements  In  the 
elastic  region  are  still  given  by  Eqs.  (4.6)  and  (4.7).  Boundary  con¬ 
ditions  for  the  problem  are 


o=0 

r 


r  =  k 


o=o 
r  r 


€r  "  €e  =  €l 


(4.18) 


r  =  b 


o  =  -p 
r  * 


(ill)  Stresses  and  Strains  In  the  Locked  Region 
From  Eqs.  (4.l6)  and  (4.17)  one  has 


—  -  e.  =  0 
dr  r  l 


(^•19) 


53ie  solution  of  this  equation  is 


u  =  C^r  +  €^r  log  r 


(4.20) 


Then, 
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er  <=  C3  +  ei(l  +  log  r) 

(*•21) 

e0  =  C3  +  £i  l0S  r 

From  Eqs.  (4.15)  and  (4.21)  one  can  write 

aQ  =  \  (3K(3C3  +  3c j  log  r  +  €£)  -  or3  (4.22) 

By  substituting  this  expression  in  the  equilibrium  equation  (4.14)  one 
obtains 

do  3o 

_E  t  .  |£  (3C3  *  3eJ  loe  r  ♦  €,)  -  0  (4.23) 


which,  when  solved,  yields  the  following  expression  for 
_  3CU 

o  =  +  3KC_  +  3Kc .  log  r 

r  3  3  l 


whence 

_  -z  Cr  €  . 

o0  -  -  I  -3  +  3KC3  +  3Kc  log  r  +  3K 

r  ' 

By  substituting  Eqs.  (4.20),  (4,24),  (4.6),  (4.7a-e), 
conditions  (4.l8b-e),  one  obtains 


(4.24) 


in  the  boundary 


C 


1 


E_  4G  }L  ll 

3K  ‘  3K  .  3  3 
b 


C 


2 


£_  4G  k3  U 
3K  "  3K  "3  3 

D 


log  k 


(4.25) 


Ck  =  +  ~  (4G  +  3K)k3e£ 
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These  constants  have  been  expressed  in  terms  of  k  and  p.  The  rela¬ 
tionship  between  k  and  p  can  be  obtained  from  Eq.  ( 4 .  l8a )  .  Thus 


.3  e.  ,3  e . 

(40  +  3K)  |  -  4G  ^  =  0 

a  b 


(4.26) 


Taking  the  limit  as  b  -»  «,  i.e.,  as  a/b,  (k/b)  -*  0  ,  Eqs.  (4.25)  and 
(4.26)  will  become 


°i  ■  -  it 


C3  =  "  3K  "  €l  108  k  "  3 

*3e. 

C]|  =  +(4G  +  3K)  -y* 


k3€/ 


(4g  +  3K)  ^  ^  +  3Ke/  log  |  =  p  +  K£f 


(4.27) 


Now,  the  expressions  for  the  stresses  are 


.3  e 

Or  =  +(4G  +  3K)  ^  -  p  +  3Ke/  log  l-K^ 

r 


oe  =  -  |  (1.G  .  3K)  ^  -  p  ♦  3^  log  §  ♦  Y 

r 


(4,28) 


Variations  of  a  ,  aa  with  the  radius  r  (for  various  values  of  k)  are 
r  w 

plotted  on  Fig.  4.2.  In  Fig.  4.3,  variations  k/a  with  p/Ke^  are 
plotted  for  different  values  of  G/K  and  b/a. 
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Fig.  k. 2b. 
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Fig.  4.3a 
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oiui 


The  expressions  (^-28)  show  that  there  is  a  stress  concentration  at 
r  »  a.  The  stress  o@  at  r  =  a  is  given  by 

°6  =  ‘  I  (p  +  3Keilog  |)  (4.29) 

When  k  a,  the  stress  concentration  factor  reduces  to  the  elastic 
value  of  3/2 . 

(iv)  Apparent  Bulk  Modulus 

One  can  measure  the  reduction  in  volume  A V  of  a  solid  sphere 
subjected  to  uniform  external  pressure  p,  by  measuring  the  reduction  in 
the  outer  radius  of  the  sphere ,  i . e . , 

AV  =-Wb2(u)  .  (4 . 30 ) 

r=b 

where  2b  is  the  diameter  of  the  sphere.  Ibis  equation  is  consistent 
with  the  assumption  of  small  displacements. 

If  the  sphere  were  made  of  elastic  material,  AV  can  be  calculated 
from  Eq.  (4.10a),  with  a  =  0,  i.e., 

AV  =  +4ub3 

Then  one  can  obtain  the  bulk  modulus  as 


On  the  other  hand,  if  the  sphere  has  a  small  cavity  of  radius  a,  and 
one  measures  the  reduction  in  volume  by  measuring  the  reduction  in  the 
outer  radius  only  and  calculates  the  bulk  modulus  from  the  preceding 
equation  one  obtains 

!.<£=&£  (4.31a) 

1  .  3a3 

*  37 
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This  value  of  the  hulk  modulus  K  Is  defined  as  the  apparent  hulk 
modulus.  K  Is  lower  than  the  true  value  of  the  hulk  modulus  K.  This 
means  that  the  presence  of  a  hole  softens  the  elastic  sphere. 

Let  us  now  consider  a  sphere  made  of  distortlonally  locking  material 
and  having  a  spherical  cavity  of  radius  a.  Then  from  Eqs.  (4.20)  and 
(4.30a)  the  value  of  the  apparent  bulk  modulus  is 


Ki  " 


1  €1  k3  ( 

*  ~7\  $ 


(4.31b) 


The  value  of  K^/K  has  been  plotted  in  Fig.  4.4  for  values  of  b/a  = 
2,3,  and  V  =  0.25.  One  can  see  that  locking  slightly  hardens  the 
material.  For  b/a  =  »,  k/b  =  0  for  finite  p  and  =  K  . 

Now,  let  us  consider  a  solid  body  of  volume  V  as  Bhpwn  in  Fig. 

4-5  with  several  small  cavities.  Let  the  volume  enclosed  by  each  cavity 
be  v,  the  number  of  cavities  n  and  the  diameter  of  the  cavity  2a.  If 
the  cavities  are  evenly  distributed  as  shown  in  the  figure,  one  can 
imagine  the  volume  V  to  be  made  of  n  fictitious  spheres  of  diameter 
2b  =  2  3V/4n.  Then  the  approximate  value  of  the  reduction  in  volume 

AY  computed  by  measuring  only  the  changes  in  the  outer  dimensions  of 
the  body  is 


AV  =  pV  [  J 

4Gb 


if  the  body  is  elastic  and 


(1  +  i| 


if  the  body  is  locked.  Then  one  can  calculate  the  apparent  bulk 
modulus  K.  The  variation  of  K  with  p  for  n  *  12,  4,  1  .  and 
V/v  =  324  are  shown  in  Fig.  4.6. 
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Fig.  4.4. 
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Fig.  4-5- 


0.5 


Fig.  4.6. 
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IV. 2.  Stresses  Around  a  Small  Spherical  Inclusion  in  a  Body  Subjected 
to  Uniform  External  Pressure 

In  this  section  we  shall  consider  a  distortionally  locking  body 
containing  a  small  spherical  rigid  inclusion  of  radius  a.  The  radius 
of  the  spherical  inclusion  is  assumed  to  be  very  small  compared  with 
the  other  dimensions  of  the  body.  If  the  body  is  subjected  to  a  uniform 
external  pressure,  the  stress  distribution  will  be  different  from  the 
uniform  hydrostatic  stress.  This  stress  distribution  will  now  be  studied. 

(i)  Elastic  Solution 

One  can  follow  a  procedure  very  similar  to  that  of  section 
IV. 1.  Then,  we  shall  again  consider  a  large  sphere  of  radius  b  con¬ 
centric  with  the  sphere  of  radius  a.  The  radius  b  is  assumed  to  be 
very  large  compared  with  a.  Often  the  difference  between  the  stress 
distribution  in  the  body  and  the  uniform  hydrostatic  stress  distribution 
will  be  very  small  at  the  outer  radius  r  =  b. 

Often,  one  can  assume  that  the  stress  distribution  at  r  =  b  »  a 
is  essentially  a  uniform  hydrostatic  stress.  After  deriving  the  formulae, 
the  expressions  are  simplified  by  considering  the  limit  as  b  tends  to 
00, 


When  the  applied  external  pressure  p  is  very  small,  the  hollow 
sphere  of  inner  radius  a  and  outer  radius  b  deforms  as  an  elastic 
body.  Then  the  general  expressions  for  stresses,  strains  and  displace¬ 
ments  are  given  by  Eqs.(4.6)  and  (4.7).  In  this  case  the  boundary 
conditions  can  be  expressed  by  the  following  equations 

r  =  a  u  =  0 

r  =  b  a  =  -p 

r 

b  ♦  » 

By  using  these  boundary  conditions,  one  can  evaluate  the  constants  C1 
and  Cg  in  Eqs.  (4.6)  and  (4.7).  They  are 


(4.32) 
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3K  +  4g  nr 
b3 


C  =  V  , 
^2  3K  '  4g 

~3  +  73 


Taking  the  limit  as  b  one  has 


c  -  -  E—  c  -  22_ 

C1  "  3K  *  C2  3K 


(4.33) 


(4.34) 


The  expressions  for  the  stresses  and  strains  are  then  given  by  the 
following  equations : 


€ 

r 


_  £_  . 

3K 


2  £- 
3K 


3 

a 

3 

r 


€ 


e 


.  £_  +  2L. 

3K  3K 


3 

a 

3 


r 


)  (4-35) 


4g  a3 

ar  =  "P  •  3K  P  “3 
r 

2G  a3 

°e  =  'p  3K  P  “3 
r 


>  (4.36) 


The  locking  condition  for  the  material  is  again  assumed  to  be  of  the 
form  of  Eq.  (4.12).  Then,  the  body  deforms  as  an  elastic  body  if 
every  where 

‘e  -  S  ■  |  <  U  '“■3T) 

r 

and  the  body  begins  to  lock  at  r  =  a  when 
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(4.38) 


£ 


K 


■  € 


l 


(ii)  Problem  After  Incipient  Locking 
If 

P>Kfe/  (4.39) 


the  body  contains  a  region  which  is  distortionally  locked.  Because  of 
spherical  symmetry  one  can  assume  that  the  region  a  <  r  <  k  is  locked 
while  the  region  r  >  k  is  still  elastic.  Hie  expressions  for  the 
stresses,  strains  and  displacements  in  the  elastic  region  are  given  by 
Eqs.  (4.6)  and  (4.7).  However,  the  stresses,  strains  and  displacements 
in  the  locked  region  must  satisfy  the  equilibrium  equation  (4. 14), 
dilatational  elasticity  (4.15),  kinematics  (4.17)  and  the  following  ideal 
locking  condition. 


€0 


€  = 
r 


(4.40) 


The  boundary  conditions  for  the  problem  are 


r  =  a 
r  =  k 


r  =  b 


(b  •*  ») 


> 


(V.W) 


(iii)  Stresses  and  Strains  in  the  Locked  Region 
From  Eqs.  (4.40)  and  (4.17),  one  has 


du  u  _  n 
d?  ■  r  +  €l  ~  0 
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(4.42) 


Hie  solution  of  this  equation  yields  the  following  result : 

\ 

u  =  C^r  -  e^r  log  r 

7r  =  C3  -  e^(l  +  log  r)  >  (4.43) 

«0  "  C3  "  ef  108  r 

Hien,  the  dilation  elasticity  Eq.  (4.3a)  yields  the  following  relationship 

between  a.  and  a  : 

y  r 

o0  =  i  [3K(3C3  -  3e£  log  r  -  e£)  -  (4.44) 

Hiis  equation  together  with  the  equilibrium  equation  (4.2)  then  results 
in  a  differential  equation  for  a^,  i.e., 

da  3a  , 

+  — -  -  (3C,  -  3e  iog  r  -  €  . )  =  0  (4.45) 

dr  r  r  3  l  t  '  " 


Hie  solution  of  this  equation  is 


_  3C4 

or  =  — j-  +  3KC3  -  3K€j  log  r 


and  \ 

3  C4  5Kel 

°e  -  -  H  +  3KC3  -  3Kez  log  r  -  -Q1 


(4.46) 


Eqs.  (4.6),  (4.7),  (4.43),  (4.46)  and  the  boundary  conditions  (4.4l) 
yield  the  following  values  for  the  constants  C^,  C^,  C3  and  C^. 


C 


1 


4G  k3  ^ 
3K  3K  75  3 

D 


L 
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C3  “  el  l0S  a 

<\-t  (-5K€/  l0s|  - » ♦  -“r)  -  -k5tsK*ko)^ 

Hie  relationship  k  and  p  is  expressed  in  the  form  of  the  following 
equation: 


3  K  3K  73  3  3 

D 


e  l  106  k 


(4.48) 


Hie  stresses  are  then 


™  £  .k 

a  =  -(3K  +  4g)  +  3K£^  log  -  a  <  r  <  k 

3r 


oQ  =  -p  +  2G 


»& 

+  3Ke^  log 

k3«j 

fe-S 

k3 

(i_  +  A 

3  €l 

\b3  r/ 

a  <  r  <  k 


k  <  r  <  b 


k  <  r  <  b 


\  (4.49) 


If  one  considers  the  limiting  values  as  b  -»  »,  the  constants  Cg, 

and  become 


C3  =  U  l0«  8 
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ck  -  r  (5KE,  1o«  I  - » - *  r) 

Then,  the  relationship  between  k  and  p  can  be  expressed  in  the  form 
of  the  following  equation 


.  i  £ 


3  K 


=  l0g  f 


The  expressions  for  the  stresses  are  then 


o  =  -  -s-j  (3K  +  4g)  +  5Ket  log  - 

r  3r  1 


a  <  r  <  k 


_  £  .k  e 

oa  =  -V-  (3K  +  4g)  +  3Ke .  log  -  -  3K  ^  a  <  r  <  k 

0  6r  £  r  2  -  - 


40  k 

ar=-p‘““€i  r>k 

r 


>  (*-50) 


x  2G  k“  . 

°e  =  "p  3  5  e/  r>k 


The  variation  of  a  and  o_  for  various  values  of  k/a  are 
r  0 

shown  in  Fig.  4.7.  b/a  has  been  assumed  to  be  ».  In  Fig.  4.8  varia¬ 
tions  of  p/Ke^  with  k/a  are  plotted  for  different  values  of  G/K 
and  b/a-  The  preceding  expressions  for  the  stresses  yield  the  following 
values  for  the  stresses  at  the  surface  of  the  rigid  inclusion. 


3 

a 


(-3K£,  l°s  I 


p 


=  -  ^  (3K  +  4g) 
3a3 


e .  e  k3  3K€ 

+  K  f-  =  -*3  (3K  +  4G)  -  -f 
6a 


If  the  body  were  completely  elastic  the  stresses  at  the  surface  of  the 
inclusion  would  be 
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,0  |  S) 


if  the  material  is  elastic,  and 


^  1  tl  *3  /■,  +  i\ 

K  ‘  P  7  \  3k/ 


(4.51b) 


if  the  material  is  locked.  The  variation  K^/K  for  different  values  of 
b/a  are  shown  in  Pig.  4. 9a.  In  Fig.  4.9b,  variation  of  K^/K  with 
k/a  for  b/a  =  3,«  are  plotted  for  the  case  of  a  spherical  cavity  and  a 
spherical  inclusion.  The  different  effects  of  the  cavity  and  the  in¬ 
clusion  can  be  seen  in  the  figure. 

IV. 3.  Stresses  Around  a  Circular  Cylindrical  Hole  in  a  Body,  in  Plane 
Strain  Condition,  Subjected  to  Uniform  Pressure  Along  the  Edges 
of  Every  Cross-Section 

We  shall  now  consider  a  body  (Fig.  4.10)  whose  length  in  the 
z-direction  is  large  compared  with  the  dimensions  in  x  and  y  directions. 
The  body  is  subjected  to  uniform  pressure  p  applied  at  the  edges  of  every 
cross-section  perpendicular  to  the  z-direction.  The  sections  normal  to 
the  z-direction  are  assumed  to  be  restrained  from  deformation  in  that 
direction.  The  stresses  in  such  a  body  with  a  circular  cylindrical  hole 
of  radius  a  will  be  investigated,  "a"  is  very  small  compared  with  the 
dimensions  of  the  body. 
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(i)  Elastic  Solution 


When  the  applied  pressure  is  small  enough,  the  body  deforms 
like  a  completely  elastic  body.  We  shall  anlayze  the  stresses  under  such 
conditions  by  considering  a  circular  cylinder  of  radius  b  concentric 
with  the  cylinder  of  radius  a.  If  b  »  a  the  stresses  at  b  are 
essentially  a  uniform  pressure  p.  Then  one  can  assume  (cylindrical 
polar  coordinates  have  been  used  for  analytical  convenience)  at 


r  =  b  a  =  -p 

r  * 

and  at 

r  =  a  a  =  0 

r 


(4-52) 


If  we  now  assume  an  axisynmetric  stress  distribution  in  the  region 

[3]  * 

a  <  r  <  b  we  have  the  general  expression  for  stresses1 


(4.53) 

Ety  using -the boundary  conditions  (4.52)  we  can  determine  the  constants 
and  Cg.  Then  the  stresses  are 


The  corresponding  strains  are 


(4.55a) 


(4.55b) 


ce 


p 


1  +  v 
E 


2.2 

a  b 

1 - 


>2  ) 


(4.55c  ) 


The  body  is  made  of  distortionally  locking  material.  In  this  case 
the  condition  of  distortional  locking  is  assumed  to  be  of  the  form 


l« 


i 


=  €l 


(4.56) 


where  e i  and 
that  \€i  -  e^l 
takes  place  at 


£.  are  to  be  chosen  out  of  three  principal  strains  such 
J 

is  a  maximum.  Fran  Eq.  (4.56)  one  can  see  that  locking 
r  =  a  when 


€r  ■  ce 


-  all  -t-  V)  p*  e 

E  b2  -  a2  1 


(ii)  Problem  After  Incipient  Locking 
If  the  pressure 


P  > 


Eei(b2  -  a2) 
2(1  +  v)  b2 


(4.57) 


(4.58) 


the  body  contains  a  region  which  is  locked  and  a  region  which  is  elastic. 
Because  of  axisymnetry  we  can  expect  that  the  region  a  <  r  <  k  is 
locked  while  the.  region  r  >  k  is  elastic.  The  expresssions  for  the 
stresses  and  strains  in  the  elastic  region  are  still  given  by  Eqs.  (4-53) 
and  (4.56)  while  the  stresses  in  the  locked  region  should  satisfy: 

(a)  The  Equilibrium  Equation 

(4.59) 


(4.60) 


do  a  -  a 

— I  +  -I - S  =  0 

dr  r 


(b)  Dilatation  Elasticity 


7r+7e  =  k  (°r+W 
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(c)  The  Ideal  Locking  Condition 


e  -  ea  =  e. 
r  9  £ 


(4.6l) 


(d)  Kinematics 


7  =  7  .  E 

r  dr  9  r 


(4.62) 


(d)  The  Incremental  Stress-Strain  Relationship  which  will 
be  derived  in  the  next  sub-section. 

The  boundary  conditions  for  the  problem  can  be  expressed  by  the 
following  equations 


r  =  k  : 


r  =  b  : 


a=0 

r 


o=o 
r  r 


cr  -  ce  =  €t 


°r  =  -p 


(4.63) 


(iii)  Incremental  Type  of  Stress-Strain  Relationship 


In  the  present  problem  the  locking  condition  is  given  by 


*  •  £e  -  er  -  cl  ■  0 


As  discussed  in  Chapter  II,  the  stress  components  o^,  oQ  and  o^  can 
be  represented  as 


a  =  a  .  +  o 
r  ri  re 


°e  =  °e/  +  °9e 


a  =  a  .  +  o 
z  zi  ze 


(4.64) 


where  o^,  and  o2^  are  the  locked  parts  of  the  stress  components, 
and  "c^,  are  the  elastic  parts  of  the  stress  components.  From 
the  Eq.  (2.27a)  we  have 


-  94  - 


&°ei  m  g?j  =  w 


Jrl 


5az/  =  0 


(4.65) 


The  relationship  between  the  elastic  components  of  stress  and  strain  are 
given  by  Hooke ' s  law ,  i .  e . , 

a  -  0o  =  2G(7  -  eQ) 

re  0e  v  r  9 


a  -a  =  -2Ge 
ze  re  r 


®«e  *  °re  +  °6e  =3K(“r  +  *9} 


From  this  equation  and  the  locking  condition  (4 . 6l)  we  have 


(4.66) 


6(7q  -  or)  =  6(7^  +  oQt  -  are  -  orJ)  =  26KT 


(4.67a) 


6a  =  5a  +  6a  .  =  6a  =  v(6a  +  6a  )  =  v(6o  +  6a  )  (4.67b) 

z  ze  zz  ze  re  We  r  a 


Equation  (4.67b)  can  be  integrated  to  yield 


°z  =  v(or  +  °e} 


(4.68) 


Another  useful  formula  can  be  obtained  by  using  (4.68)  and  the  dilata- 
tional  elasticity  Eq.  (4.60) 


;r  +  °e =  rrr  (rr  +  re> 


(4.69) 


(iv)  Solution  in  the  Locked  Region 

From  Eqs.  (4.6l)  and  (4.62)  we  get  the  following  differential 
equation  in  u. 
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(4.70) 


du  u  ^ 

dr  "  r  '  ei  =  0 


which  when  solved  yields  the  following  result. 


u  =  C^r  +  c^r  log  r 
er  =  C3  +  <l l  (1  +  log  r) 

"e  =  C3  +  el  log  r 


Then  from  Eqs.  (4.69)  and  (4.71)  we  have 


(4.71) 


1  +  v  ^2C3  +  2ei  lo8  r  +  e^)  - 


(4-72) 


The  equilibrium  equation  (4.59)  then  reduces  to  a  differential  equation 


in  o 


do  2o 


IT  +  IT  =  1TTT} 


(4.73) 


Hence , 


3K  „  .  3K  ,  4 

°r  =  TTU  c3  +  TT~v  €i  log  r  +  “5 

r 


(4.74) 


(v)  Matching  the  Boundary  Conditions 

The  solution  in  the  elastic  region  is  given  by  Eqs.  (4.53)  and 
(4.56).  By  using  the  boundary  conditions  (4.63)  we  find 


C  -  ■  E  k2e 

1  2(1  +  V)  £ 


(4.76a) 


2Ti~TvJ  ^  €l  '  p 


(4.76b) 


v  1^2v  k2 

-*t  log  k  -  | - ^ - ■£  e£ 

~  b 


-  -^loga-e^l-v)^ 
a 

(4.76c) 


BeJ(l-  v) 

C4  *  (1+  v)(l-  2V ) 


(4.76d) 


and  the  relationship  between  k  and  p  is  given  by  the  following 
equation 


P  = 


TT 


26 J  f/,  ..s  k2  a  1  1  -  2V  k2! 

7  v)(l  -  2v)  [_(!  -  v)  7  +  108  k  ■  2  ’  “ - 


(4.77) 


Now  if  we  consider  the  limit  as  b  -*  <»  the  expressions  for  C^,  C 2>  C3 


and  become 


C,  = 


1  =  2(l  +  v) 


k2e, 


c2  =  -P 


Be/1  -  v)  2 

C4  =  (l  +  v)(l  -TvJ  k 


Then  the  stresses  are  given  by  the  following  equations : 


k2 

°r  =  2(l  +  v)  2  "p 


^  t  k2 

ae  =  '  2(i  +  v;  r2  _p 


r  =  (l  +  V)(l  -  2vJ  [j1  ’  2  +  log  k  “  l]  "p 

;e  =  -  (i  +  v)(i  :"2vl  [(1  •  v)  7  "  log  f  '  i]  _I 


(4.78) 


where  k  follows  from 
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p  = 


Be 


(1  +  v)(l  -  2v) 


-  V)  %  +  log  & 
a 


The  stress  aQ  at  r  =  a  is 


K £.  . 

-2*  +  irrv 7  log  ! 


(4-79) 


This  is  higher  than  the  value  of  o@  derived  on  the  assumption  of  elastic 

behavior  of  the  material.  In  Fig.  4.11,  the  variation  of  a  and 

r0 

with  r  Is  shown  for  different  values  of  k.  b/a  has  been  assumed  to  be  00 . 
In  Rig.  4.12a,  p/E  has  been  plotted  against  for  different  values  of  v  . 

The  effect  of  different  b/a  is  shown  in  Fig.  4.12b.  In  concluding 
this  section,  one  can  observe  that  the  problem  of  computing  stresses 
around  a  circular  cylindrical  rigid  conclusion  in  a  body,  in  plane 
strain  conditions  subjected  to  uniform  pressure  along  the  edges  of  every 
cross-section  can  be  worked  out  in  a  similar  way. 

IV. 4.  An  Example  to  Illustrate  the  Limitations  on  the  Displacement 
Boundary  Condition 

It  has  been  proved  in  section  1-13  that  if  the  displacements  are 
prescribed  on  the  surface  of  the  body  so  that  they  increase  from  zero 
value  proportionally,  the  prescribed  displacement  can  be  increased  cnly 
up  to  that  value  which  makes  the  body  completely  locked.  This  will  now 
be  illustrated  by  a  simple  example. 

We  consider  a  hollow  sphere  of  inner  radius  a  and  outer  radius 
b  with  prescribed  displacement  u  =  xCi  at  r  =  a  and  u  =  n0  at 
r  =  b.  a.  is  assumed  to  be  greater  than  The  prescribed  displace¬ 
ments  are  assumed  to  increase  from  their  zero  value  proportionally,  i.e., 
n  increases  from  its  value  zero  to  its  present  value.  For  small  values 
of  n  the  sphere  deforms  like  a  completely  elastic  body.  Further,  this 
problem  has  spherical  symmetry.  Then  the  stresses,  strains  and  dis¬ 
placements  are  given  by  Eqs.  (4.7)  and  (4.6).  By  using  the  boundary 
conditions  the  constants  and  Cg  can  be  evaluated.  Thus, 
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(Ik  80) 


„  2,2  Ob  -  Pa 

C2  «  Mb  -y-  ^ 

b  -  a 


_  _  na  nb  Oft)  -  Pa 

C1  “  a  "  a  73  3 

b  -  a 


J 


The  stresses,  strains  and  displacements  in  the  body  can  then  be  calculated 
using  these  values  of  and  Cg.  The  sphere  is  capable  of  locking  dis- 
tortionally.  The  locking  condition  is  assumed  to  be  of  the  form 


<€e 


(4.81) 


From  Eqs.  (4.7a)  and  (4.87)  -  €  can  be  computed.  Its  value  is 

w  r 


,  3,2  a  a  ‘  P 

€e  "  r  =  3na^  73 - r 

d  -  a 


1 

“5  ’ 

r 


a  >  p 


(4.82) 


This  has  the  maximum  at  r  =  a.  Then  the  sphere  begins  to  look  at 
r  =  a  when 


n 


(4  83) 


(i)  Problem  After  Incipient  Locking 

As  n  is  further  increased  from  its  value  given  by  (4.83)  the 
sphere  will  contain  a  region  which  is  locked  and  a  region  which  is 
elastic.  From  spherical  symmetry  we  can  assume  that  the  region 
a  <  r  <  k  is  locked  and  the  region  k  <  r  <  b  is  still  elastic. 

In  the  region  a  <  r  <  k  the  general  expressions  for  stresses, 
strains  and  displacements  are  given  by  Eqs.  (4.20),  (4.21),  and  (4.24) 
with  replaced  by  i.e., 
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u  *  C3r  -  c^r  log  r 

«r  =  °3  '  €i(l  +  106  r) 
7e  a  c3  '  log  r 


_  3C^ 

o  =*  +  3KC_  -  3Ke .  log  r 

r  3  3  l 

r  J 


(4.84) 


Hie  general  expressions  for  stresses  in  elastic  region  are  again  given 
by  Eqs.  (4.6)  .  Now,  the  boundary  conditions  can  be  expressed  by  the 
following  equations 

r  =  a  u  =  net  1 


r  =  b 
r  =  k 


u  =  ng 
u  =  u 


a 

r 


=  a 

r 


e 


r 


=  -€  , 


Then 


(4.85) 


°2  •  *1  r 

C3  ■  f  *  el  l0«  8 


(4.86) 


(40  -  3K) 


Hie  relationship  between  n  and  k  is  given  by 


f 


log  a 


€  ^  log  k 


(4.87) 


102  - 


The  sphere  will  be  completely  locked  when  k  =  b,  i.e.,  when 


n 


b 


e£ 


(4.88) 


(ii )  Problem  After  the  Sphere  is  Completely  Locked 

Let  us  now  attempt  to  study  the  sphere  in  the  completely 
locked  state.  In  this  state  the  expressions  for  stresses,  strains  and 
displacement  throughout  the  sphere  are  given  by  Eqs.  (4.84).  The  boundary 
conditions  are 

r  =  a:  u  =  nG 

r  =  b:  u  =  n0 

From  (4.89a) 

C3  =  log  a  +  r 

From  (4.89b) 

c3  =  et  lo« b  + 

Contradicting  conditions  (4-90)  and  (4-91)  for  b  require;  that 

€l  log  t  "  n  (b  - 

This  can  be  reconciled  for  any  n  only  if  a  =  b.  Therefore,  the  pro¬ 
portional  increase  of  surface  displacement  is  not  possible  for 


(4.89a) 

(4.89b) 

(4.90) 

(4-91) 


n  > 


v  b 

£lb  lQg  a 
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V.  WAVE  PROPAGATION  IN  LOCKING  MEDIA 


V.l  One -Dimensional  Wave  Propagation 

(l)  Stress -Strain  Relationship  In  the  Uniaxial  Stress 
System  and  Volumetric  Locking 

In  this  section  we  shall  study  the  wave  propagation  under  a 
uniaxial  stress  system.  As  a  first  step  toward  the  study,  the  stress- 
strain  relationship  in  a  volumetrically  locked  region  applicable  in 
this  case  will  be  derived. 

The  locking  condition  for  the  material  is  assumed  to  be 
€+€+£=£.  (5-1) 

x  y  z  1 


After  locking,  the  deviatoric  parts  of  stress  and  strain  tensors  are 

still  related  by  Hooke's  law.  In  the  one-dimensional  stress  system  the 

stress  tensor  consists  of  only  one  component  ox  .  The  strain  tensor 

contains  three  components  of  normal  strain  e  ,  e  ,  €  .  If  the  devia- 

x  y  z 

toric  parts  of  strain  and  stress  tensors  are  obtained  and  Hooke ' s  law 
is  used  one  obtains 


2a  =  2G(2l-  e  -  e  ) 
x  x  y  z 

-a  =  2G(-I  +  2e  -  e  ) 
x  x  y  z 

-a  =  2G(-e  -  e  +  2e  ) 
x  '  x  y  z' 


(5.2) 


Only  two  of  these  three  relations  are  independent.  By  simplifying 
these  equations  one  obtains 


e 

y 


e 

z 


(5.3) 


Eqs. (5.3)  together  with  the  locking  condition  yield  the : following 
relationships 


3Ge  -  Gc . 
x  l 

H  ax 

3  "  5g 


(5.4a) 

(5.4b) 
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These  equations  express  the  relationship  between  ax,  e*  and  e  in 
the  locked  region. 

The  locking  condition  (5*1)  can  be  expressed  in  terms  of  stresses 
also.  If  the  material  is  elastic 


1  -  2v 

e  +  £  +  e  =  — - —  (o  +  o  +  a  ) 

x  y  z  E  '  x  y  z' 


In  this  case  o  =  a  =  0 

y  z 


Then  Eq.(5.l)  becomes 
E€. 


ax  ~  1  -  2v 


(5.5) 


Then  one  can  state  the  following  inequalities.  The  material  is  elastic 
if 


Ee, 


a  ^  - — 

x  1-  2v 


and  locked  if 


Ee, 


ax  ^  1  -  2v 


The  stress-strain  diagram  is  as  shown  in  Fig. 5.1 


Fig.  5-1. 
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The  straight  lines  OA  and  AB  are  defined  by  the  following  equations 


Let  us  consider  a  bar  of  uniform  cross  section  made  of  elastic¬ 
locking  material.  It  is  assumed  that  the  average  stresses  on  any  cross 
section  can  be  approximated  by  a  one -dimensional  stress  system.  Then, 
one  can  apply  Newton's  law  of  motion  to  an  element  of  bar  of  length  dx 
as  shown  in  Fig. 5*2  .  The  displacement  in  the  x-direction  is  represented 
by  u  and  the  velocity  by  v  =  du/dt  . 


v  v  +  dv 


Fig.  5-2. 


In  further  discussion  of  one-dimensional  wave  propagation  the 
usual  notation  of  barred  letters  to  denote  the  stresses,  strains  and 
displacements  in  the  locked  region  will  not  be  applied  because  we  are 
working  most  of  the  time  only  with  the  stresses  in  the  locked  region. 
The  equation  of  motion  can  then  be  written  as 


*ax  dv 


*57“  p  Jt 


In  general  the  stress-strain  law  can  be  expressed  as 
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0  -  f  (c) 


(5.6) 


If  one  defines  the  normal  force  N  such  that 
N  =  A  a 

x 

and  the  mass  per  unit  length  4  as 
4  *  Ap 

the  equation  of  motion  and  the  stress-strain  relationship  can  be  re¬ 
written  in  the  following  way 

N'  =  uu  ,  N  =  Af (u ' )  (5-7) 

where  (  ) 1  represents  partial  differentiation  with  respect  to  r  and 
(  )  represents  partial  differentiation  with  respect  to  t  .  Elimina¬ 
tion  of  N  from  Eq. (5-7)  yields 

a£tu"  =  4U  (5.8) 


By  introducing  the  symbol 


E(u') 


df 

du' 


equation  (5-8)  becomes 


u 


If 


AE(u')  U 


(5.9) 


From  this  equation  one  can  infer  that  there  is  a  wave  of  velocity  c 
such  that 


c 


2  _  AE(u') 
4 


(5.9a) 


moving  in  the  locked  region. 
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Alternately,  one  can  derive  the  equation  of  motion  by 
the  discontinuities  at  the  wave  front. 


Fig.  5-3. 

By  the  impulse  momentum  theorem  one  can  write  (See  Fig. 


dtA(°xl-ax2)  =  pA^2"  V3* 


i.e., 

ax2 "  Vpc(V  V 


or 

N1  "  N2  =  pc^2"  “l^ 


where 


By  introducing  the  Jump  symbol 
A(  )  =  (  )g  -  (  \ 

one  can  rewrite  equation  (5.10a)  in  the  following  way 
UcAu  *  -  AN 


considering 


•3) 


(5.10) 


(5.10a) 


(5.11) 
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From  the  condition  of  continuity  of  displacements  along  the  wave  front, 
one  can  write 

Ug  -  =  Au  =  0 

By  differentiating  this  equation  along  the  wave  front  in  the  x-direction 
one  obtains 

S  (Au)  =  [h  +  c  Au  =  Au'  +  c  =  0  (5.12) 

Eqs.  (5.11)  and  (5*12)  can  be  simplified  to  obtain 

4C2Au  ’  »  AN  (5.13) 

If  one  makes  the  assumption  that  the  discontinuities  Au,  AN,  Au' 
carried  by  the  wave  are  infinitesimal  quantities,  one  can  write 

AN  =  AAf(u')  =  AE(u ' )Au ' 

Then, 

UC2Au'  =  AE(u')Au' 


i  •  6  •  y 


The  result  represented  by  Eq.  (5.1*0  is  the  same  as  that  obtained  by  the 
differential  equation.  However,  by  the  study  at  wave  front  we  have 
learned  that  the  results  (5.1*0  and  (5* 9a)  are  true  only  for  infinites¬ 
imal  discontinuities.  Then  the  next  problem  will  be  to  study  whether 
there  is  a  possibility  of  the  discontinuity  being  finite.  As  a  first 
step  toward  the  study  the  variation  of  discontinuity  during  the  propa¬ 
gation  of  the  wave  will  be  studied  in  the  next  section. 
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(ill)  Variation  of  Discontinuity  During  the 
Propagation  of  the  Wave 

By  differentiating  Eq. (5.11)  along  the  line  of  propagation 


of  the  wave  in  the  x-direction,  one 

nc  +  nc  'Au  +  ^  Au  =  -AN'  -  i  Afl  (5.15) 

From  Eq.  (5*7)  one  can  write  that 

AN '  =  HAU  (5.16a) 

A&  =  A(EAu ' )  (5.1.6b) 

Eq. (5.16b)  can  then  be  written  as 

AN  =  A  u'Au 1  +  pc2Au'  (5.17) 


This  equation  is  for  infinitesimal  waves  only.  Further,  from  Eq. (5.9a) 
one  has 


2cc ' 


A  _dE 
U  du' 


u 


2cc 


A  _dE 
U  du' 


(5.18) 


By  substituting  (5.16a),  (5. IT) ,  and  (5.18)  inEq.(5.15),  the  following 
equation  is  obtained: 


^  s  wo  ♦  b  3- 1»"  *  J  » *  £■  4'to'  -  0 


2c  du 


c  du' 


Then,  by  using  Eq. (5*12)  this  equation  can  be  simplified  to  the  follow¬ 
ing  form 


D_ 

Dx 


(Au) 


A  dE  /  ft  1  •  |\a< 

+  ~ 2  d^  (u  '  “  u  )Ab  *  0 


Upc 


(5.19) 
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This  is  a  differential  equation  in  Au  .  It  will  be  used  in  the  next 
sub-section  to  study  the  effect  of  variation  of  discontinuity  as  the 
wave  propagates  into  the  bar.  It  should  be  noted  that  this  equation  is 
valid  only  for  infinitesimal  waves. 

(iv)  Case  of  Unloading  Wave 

Let  us  consider  a  bar  prestressed  to  N  =  P  .  Let  a  stress - 
discontinuity  of  the  magnitude 

N(0,t)  =  -  Pl(t) 

be  sent  through  the  bar  from  the  end  x  =  0  .  In  this  equation  l(t) 
represents  the  unit  step  function,  i.e. ,  the  force  P  is  suddenly 
released  at  the  end  x  =  0  .  Further,  it  is  assumed  that  the  discon¬ 
tinuity  -  P  is  made  up  of  elementary  waves  of  magnitude  -  AP  which 
propagate  successively  through  the  bar. 

At  first  u(x,t),  u(x,t)  are  equal  to  zero  throughout  the  bar. 

Then,  the  second  term  in  Eq.(5.19)  is  zero.  Therefore,  Au  is  constant 
throughout  the  bar;  i.e.,  as  the  first  elementary  wave  propagates  through 
the  bar,  the  wave  creates  a  state  of  constant  velocity  and  constant  stress 
behind  the  wave.  Then,  one  can  conclude  that  u"  and  d'  are  zero 

behind  the  wave.  Thus,  for  the  second  elementary  wave  of  magnitude  -AP  , 

the  second  term  of  the  Eq.  (5-19)  is  again  zero.  The  Eq.  (5.19)  again  yields 
the  same  result,  i.e.,  Au  *  constant  throughout  the  bar.  That  is,  the 
propagation  of  the  second  elementary  wave  through  the  bar  also  creates  a 
state  of  constant  stress  and  constant  velocity  behind  the  wave. 

Thus,  one  can  see  that  all  elementary  waves  spread  through  the  bar, 

similar  to  the  linearly  elastic  waves.  However,  they  propagate  with 

different  speeds  as  given  by  Eq. (5 -9s).  It  means  that  after  the  nth 
elementary  wave  has  passed  through  an  arbitrary  point  x^  of  the  bar, 
the  stress  at  the  point  is  P  -  nAP  .  The  speed  of  propagation  of  the 
nth  elementary  wave  is 


where  Er  is  the  slope  of  the  stress -strain  curve  evaluated  at 
rAP 

a  ■  o  =  — — 
x  xn  A 


i  •  e  •  , 


From  Fig.  5.1 


Ee, 


En“3G  lf  °x  > 


Ee 

E  lf  °x<Trfe 


Also,  3G  >  E  .  Then,  one  can  see  that  the  unloading  wave  falls  apart 
into  two  waves  with  speeds  |^3g/p‘  and  j/e/p  . 

Further,  according  to  Eq.(5«H)  the  elementary  waves  add  a  con¬ 
tribution 


Au  » 


AN 

\ic 


to  the  velocity  u  .  Therefore,  the  particle  velocity  at  a  point 
allowing  an  axial  force  P  -  nAP  is 


rAP 


(5.21) 


From  Eq. (5.21)  it  is  evident  that  u(x,t)  is  a  function  of  the  stress 
ox  only.  Thus  it  does  not  depend  on  the  stress  distribution  along  the 
length  of  the  bar. 

If  one  considers  the  upper  limit  of  the  integral  in  Eq. (5*21)  as 


nAP  -  P 
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one  obtains  the  velocity  at  the  end  of  the  bar  immediately  after  the 
application  of  the  discontinuity. 

(v)  Case  of  the  Loading  Wave 

Let  us  now  consider  the  unloaded  bar  x  >  0  .  Let  a  stress 
discontinuity  of  intensity 

N  =  Pl(t) 


be  sent  through  the  bar  from  the  end  x  =  0  .  Again,  one  can  consider 
the  stress  discontinuity  to  be  made  of  elementary  waves  of  intensity 
AP  .  The  first  of  these  waves  travels  with  a  speed  c  =  c(0)  while 
the  final  elementary  wave  travels  with  a  speed 

c  =  c(P/A) 

(see  Fig. 5 -la),  which  is  greater  than  that  of  the  first  wave,  i.e.,  the 
last  wave  catches  up  with  the  first  at  a  certain  point.  Similarly,  all 
the  elementary  waves  introduced  at  a  later  time  can  catch  up  with  those 
introduced  when  nAP  <  Ee^/l-2v  .  However,  when  the  two  waves  catch  up 
with  each  other  they  cannot  pass  one  another  as  two  waves  approaching 
from  opposite  directions,  in  an  elastic  bar.  In  this  case  the  wave  that 
overtakes  passes  from  a  domain  of  higher  stress  to  a  lower  stress. 

Hence  it  becomes  the  slowest  wave.  Therefore,  one  can  see  that  the  waves 
do  not  separate  after  they  catch  up  with  each  other.  They  move  together 
with  a  common  speed,  i.e.,  the  discontinuity  at  the  wave  front  starts 
accumulating  and  becomes  finite.  One  can  no  longer  apply  the  formulae 
derived  on  the  assumption  of  an  infinitesimal  discontinuity.  In  order 
to  find  the  speed  of  propagation  of  the  wave  carrying  a  finite  discon¬ 
tinuity,  one  must  use  Eq.(5»13)  which  is  applicable  for  finite 
discontinuities,  i.e., 


2 

c 


VN2 

H(u^-  Ug ) 


1 

P 


(5-22) 


In  this  equation  and  in  further  analysis,  the  subscript  x  of  ox  and 
cx  has  been  dropped. 
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In  concluding,  one  can  see  that  the  case  of  loading  produces  a 
phenomenon  which  is  quite  different  from  that  of  unloading.  Stress  and 
velocity  discontinuities  in  the  case  of  loading  do  not  fall  apart  as  in 
the  case  of  unloading.  They  move  together  as  a  finite  wave  (called  the 
shock  wave).  The  speed  of  this  wave  at  any  point  is  given  by  Eq.(5.22). 

(vi)  Impossibility  of  an  Unloading  Shock 

The  relationships  (5*13)  and  (5.22)  are  valid  whether  the 
load  increases  or  decreases.  The  concept  of  loading  shock  has  been 
explained  on  the  basis  of  these  equations.  Then,  one  may  ask  the 
question  whether  an  unloading  shock  exists?  If  not,  what  are  the  reasons 
preventing  the  existence  of  an  unloading  shock?  This  can  be  explained 
by  energy  considerations. 

Let  us  consider  the  loading  wave.  When  the  stress  changes  from 
C1  to  °2  and  the  velocity  from  u^  to  fig  ,  the  change  in  kinetic 
energy  is 

1  .  ( *2  .2  \ 

Ux  =  g  udx  I  Ug-  u1  1 

The  increase  in  elastic  energy  during  the  same  process  is 


Ug  =  dx 


N(u')du' 


Net  work  done  during  the  process  is 
U3  =  -  (N^-  NgUg)dt 

Then  one  can  define  the  energy  loss  U  as 

U  =  increase  in  kinetic  energy  +  increase  in  elastic 
energy  -  work  done 


i.e. , 


Udx 


2 

(v  +  N(u')du’  +  ^(N^-  NgUg)dx 


(5.23) 
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The  o-€  curve  is  of  the  shape  shown  in  Fig. 5.4 


a 


Fig.  5.4. 

When  the  curvature  is  concave  upwards  the  integral  can  be  estimated 
in  the  following  way: 


A 


|  (a^  a2>  (e2~  €1>  ‘  C 


(5.24) 


where  C  is  the  shaded  area  in  the  figure.  Further,  C  is  a  positive 
quantity. 

Also,  from  Eq.(5.22)  one  can  write 


2 

pc 


t)A 


where 


n  = 


i.e. 


> 


e 


2 
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Therefore , 


) 


C 


Further,  from  the  figure  one  can  write 


°1  =  €0) 

a2  =  £0} 


Then,  the  energy  loss  U  is 


(5.25) 


The  terms  in  the  square  brackets  cancel  each  other  because  of  continuity 
of  displacements .  Then 


U  =  -  C  (5.26) 

This  means  that  there  is  more  energy  than  needed  which  makes  the  process 
possible. 

For  an  unloading  wave  the  signs  of  all  energy  terms  are  reversed, 
including  that  of  the  contribution  C  derived  from  the  nonlinearity  of 
the  stress-strain  relationship.  Instead  of  a  lose  we  would  then  have 
mechanical  energy  produced  from  nothing.  This  shows  that  upon  unloading 
the  wave  front  disperses;  a  shock  is  not  possible. 

(vii)  Interaction  of  two  Waves  Approaching 
from  Opposite  Directions 

Let  us  consider  two  waves  A  and  B  approaching  from 
opposite  directions.  The  vave  A  is  carrying  a  stress  discontinuity 
and  velocity  discontinuity  v^  .  The  wave  B  is  carrying  a  stress 
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Further,  it  is 


discontinuity  <jg  and  the  velocity  discontinuity  Vg  . 
assumed  that 


°1* 


°2  *  1  -  2v 


(5.27) 


but 


°1  +  °2  >  1  -  2V  £2 

Then  the  waves  A  and  B  are  elastic  waves  moving  with  the  elastic 
wave  speed  c  =  Je/ p  .  Further,  it  is  assumed  that  the  region  through 
which  neither  of  the  waves  has  passed  is  stress  free.  Then  from  Eq.(5.10) 
one  can  write  the  following  relationships  between  stress  and  velocity 
discontinuities 


ffl  =  pcvl 


a2  =  pcv2 


(5.28) 
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XI  X 

Fig.  5-7- 
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The  x  - 1  diagram  for  the  propagation  of  the  wave  has  been  indi¬ 
cated  in  Fig. 5 *7*  At  a  certain  point  x  ■  x1,  t  *  ^  the  two  waves 
meet.  Then  two  different  waves  D  and  E  moving  with  speeds  c^ 
and  c^  emerge  from  the  point  as  indicated  in  Figs. 5.6  and  5.7.  The 
stress  and  velocity  discontinuities  carried  by  these  waves  can  be  cal¬ 
culated  in  the  following  way; 

By  applying  Eq. (5*10)  to  the  waves  D  and  E  one  obtains 


°3  ■  ai  =  *  pc3(V  Va) 

°3  •  °2  =  pC4(V  V3} 

By  using  Eqs.(5.28),  these  equations  can  be  written  as 


(5.29a) 


°3  ‘  pcvl  =  -  PC3(V1+  V3) 

°3  •  pcv2  =  pC4^V2~  V3> 
Further,  from  Eq. (5.22) 

2  1  °1  ‘ °3 

C  SB  —  -  ■  ■ 

3  P  -  e3 

2  1  °2  ~°3 

C4  =  p  e2-€3 

But  from  Eq.  (5. 4)  and  Hooke's  law 


(5.29b) 
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Then 


2  1  al  ‘  ff3 

c_  =  — 


'3  p  ox  a3 


3G  3 


2  1  °2  ‘  °3 

c4 


P  °o 


_3  _£ 

3G  "  3 


(5,30) 


Equations  (5*29)  and.  (5-30)  are  four  equations  for  four  unknowns,  v^, 
c3  and  c^. 

In  the  special  case  when  o1  =  Og  and  =  Vg,  the  right-hand 
sides  of  Eq.  (5-30)  become  equal,  indicating  that  c^  =  -  c^.  Hie  minus 
sign  has  been  chosen  because  the  waves  D  and  E  are  moving  in  opposite 
directions.  Equations  (5-29)  can  then  be  written  in  the  following  way: 


°3  "  pCVl  *  •  pC3(vl  +  V3) 
a3  *  pCVl  =  pC4(vl  •  V3} 


whence 


v3  =  0 


-c,  =  +  c,  ,= - c 

3  4  pv^^ 


Hien,  from  Eqs.  (5*30)  and  (5 -29a) 
2  °5<V1*V3) 


’3  !i  !i  li 

E  ‘  3G  ‘  3 


(5.30a) 


By  using  Eqs.  (5.30a)  this  equation  becomes 
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(5.30b) 


'3  a 


1  p  ,  v  ei 

v,  "  K  C"  C3^“  3v, 


However,  from  Eq. (5.28) 


Then,  Eq. (5.30b)  becomes 


C3  +  C3  l“2  C  -  C 
'c 


.±) 


where 


_  /3G 
v  p 


Is  the  wave  velocity  Inside  the  locked  region  and  c  =/e/p'  is  the 
elastic  wave  velocity.  If  c^  is  determined  by  solving  the  quadratic 
equation,  can  be  obtained  from  Eq. (5- 30a) .  However,  the  quadratic 

equation  has  two  roots.  We  should  investigate  to  see  which  of  the  two 
roots  is  the  correct  answer. 


c3  can  be  expressed  as 

a  /?  -2 

C3  =  2  *  V  T  +  C 


(5.30c) 


where 


e  ,c2  -2 

l  c 

c  +  - g  C 

1  c 


30  .  £ 
p  E 


3 

2(l+v) 


121  - 


Also,  we  have 


Eei 

al  *  1  -  2v 


i.e. , 


Ee , 


3pv1  S  3pv1(l-2v) 


i.e. , 


€i  *  1-  2v 
3Vl  3c 


or 

61C  >  1  -  2v  5 
3Vl  2  3  ‘  2 ( 1+V )  C 

Then  the  minimum  value  of  a  is 

A  .  1  -  2v  3  \ 

“min  "  C  V1  +  2(l+v)  '  2 ( 1+V )  J  =  ° 

Therefore 

a  i  0  . 

The  equal  sign  is  applicable  where  =  Ee^/l-2v  i.e.,  in  incipient 
locking.  Therefore,  the  plus  sign  of  the  radical  in  Eq. (5- 30c)  gives  a 
value  of  c3  higher  than  c  .  This  is  not  possible  because  c  corre¬ 
sponds  to  the  slope  of  the  straight  line  AB  in  Fig. 5.1  or  Fig. 5*^, 
while  Cj  is  the  slope  of  one  of  the  secants  such  as  IE  or  OC  . 

Then  the  correct  c^  corresponds  to  the  minus  sign  of  the  radical  in 
Eq. (5.30c) . 
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If  in  particular 


01  =  1-  2v 

Cj  becomes  equal  to  ySG/p  because  a  =  0  .  Then 
c3  =  v1(cp  +  y/3Gp)  • 

(vii)  Interaction  of  Two  Waves  Moving  in  the  Same  Direction 

Nov,  let  us  consider  two  waves  A  and  B  as  shown  in 
Fig. 5* 8a  and  moving  in  the  same  direction.  The  wave  B  propagates 
into  the  undisturbed,  stress  free  bar.  It  carries  a  stress  discontinuity 

Es, 

°1  =  ~ 2v 

and  a  velocity  discontinuity  . 

The  wave  A  is  assumed  to  carry  a  stress  discontinuity  of  the 
magnitude  a g-  and  to  move  with  a  velocity  c^  .  The  value  of 

is  assumed  to  be  greater  than  .  Then  one  can  write  the  following 
relationships  between  stresses  and  velocities  (see  Eq. (5-10)) 

al  =  pcvl 

a2  -  01  =  pc^Vg-  v2) 

where  c  is  the  elastic  wave  speed  and  c^  can  be  calculated  from  the 
formula  (5.22). 
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Fig.  5-8b. 
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The  x,t  diagram  is  as  shown  in  Fig. 5* 8c*  The  wave  A  meets  the  wave 
B  at  a  certain  point  x  =  x^,  t  »  t^  .  Then  two  waves  D  and  E 
emerge  from  the  point.  The  stress  discontinuities  a3  -  o2  and  o3 
carried  by  these  waves  will  now  be  calculated. 

By  applying  Eq.(5.10)  at  the  wave  fronts  D  and  E  ,  one  can 

write 


°3  *  pc3v3 
°3  ”  °2  " 


(5.31) 
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Further,  from  the  equation  (5-22) 


2  1 

°3 

C-  s  “ 

3  p 

2  1 

a3  ^  ei 

3G  3 

°2  °3  3G  2  -2 

(5.32) 

C4  =  p 

(a2  -0,)  1_  p  C1 

d  5  3G 

Eqs. (5*31)  and  (5-32)  can  be  rewritten  as 

°3  = 

pC3V3 

(5.33a) 

°3  ~°2  = 

(v2-  v3)  j3Gp 

(5.33b) 

2 

C3  = 

1  °33G 

p  °3+  06  / 

(5.33c) 

Eqs. (5-33)  are  three  equations  for  three  unknowns  and  , 

which  can  then  be  determined.  Then  we  know  the  stress  discontinuity, 
the  velocity  discontinuity  and  the  speed  of  the  wave  after  intersection. 

From  Eqs. (5.33a)  and  (5«33c) 


2 

C3 


o?  .  o*3G 

-T2  -  p  ^7 


°3  +G€i 
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Prom  Eq. (5.33b)  and  this  equation 

2  .  2 

°3  +  °3 Gef  "  (®3 "  °2  '  ^Gp  v2) 

i  •  0  •  j 

+3GpVg  +2v2o^3Gp 

Q  “  ~ 

Gci  +  2a2  + 2vg /3Gp 

The  right-hand  side  of  this  equation  contains  known  quantities.  Hence, 
Oj  can  be  calculated.  Then  c^  and  can  be  obtained  from 
Eqs. (5.33c)  and  (5* 33a). 


v.2  Spherical  Wave  Propagation  under  Volumetric  Locking 

(i)  Equations  of  Motion  and  the  Speed  of  the  Wave 

In  this  section  the  dynamic  response  problem  of  an  infinite, 
homogeneous,  isotropic  and  volumetrically  locking  medium  will  be 
studied  under  spherically  symmetric  conditions.  One  can  write  the 
following  equations  governing  the  problem.  In  these  equations 
or,  ir,  u  are  the  normal  stress,  normal  strain  and  displacement  in 
the  radial  direction  r  .  Og,  are  the  normal  stress  and  normal 
strain  in  any  direction  in  the  plane  perpendicular  to  r-direction. 

The  following  equations  govern  the  motion: 

a)  Newton ' s  Law 


(5.34) 


b)  The  Stress-Strain  Law:  The  material  is  assumed  to  be 
capable  of  non-ideal  locking.  The  results  for  ideal  locking  will  be 
derived  as  a  limiting  case.  Prom  Eq. (2.36)  and  distortional  elasticity 
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'°r+2'a6  -  3K€j  +  \(€r  +  2€fl.ei)  ,  ir  +  260  >6i 
5r  -  =  20(Ir-  lQ) 


(5.35a) 


If  one  introduces  a  potential  if  such  that 


Eq. (5.37)  can  be  rewritten  in  the  following  way 


(5-38) 
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where 


2  X.  +UG 
C1  "  ~ 


and 


is  the  Laplacean  operator  under  spherically  symmetric  conditions. 

Now,  one  can  see  that  Eq.(5.38)  represents  a  wave  equation  with 
the  wave  speed 


c 


1 


(5.39) 


One  must  remember  that  the  speed  of  the  wave  is  given  by  this  equation 
only  if  the  wave  is  propagating  in  a  medium  which  is  already  locked. 

If  the  material  is  an  ideally  locking  material,  \  becomes  infinity. 
Then  the  speed  of  the  wave  becomes  infinity. 

Alternately,  one  can  derive  the  equations  governing  the  problem  by 
considering  the  discontinuities  at  the  wave  front.-  By  applying  the 
impulse -momentum  theorem  one  can  arrive  at  the  following  relationship 
(Fig. 5.9): 

A 
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(in  this  figure  AOB  is  a  cone  with  vertex  at  the  origin,  and  vertex 
angle  being  equal  to  d 0  .  The  wave  CD  sweeps  through  the  radial 
distance  r2~  ri~  dr  in  time  •) 

2  2 
T  +  T 

(N2-  Nx)dt  =  pir  —  2  2  *2)  (5.&0) 


where  and  are  the  resultant  forces  and  u^  and  ug  are  the 
radial  velocities  on  sides  (l)  and  (2)  of  the  element  shown  in  Fig. 5*9* 
Further, 


arl 

°r2 


Trr2d0 

Trr|d6 


2 

2 


Then  Eq. (5.^0)  can  be  rewritten  as 


(°r2~  °rl)  r2d02dt  +  or2Tr(2r1drd02 + dr2d02)dt 
=  (ux-  u2)pirr2d02dr(u1-  ug) 

+  pTTd02dr2(r1+ u2) 

o  2 

By  dividing  this  equation  by  d0  dt  and  by  considering  the  limit  as 
dr  and  d0  tend  to  zero,  this  equation  becomes 

ar2  “  °rl  =  pc^l'  V  (5-M) 


where 


c 


dr 

dt 


Further,  the  displacements  are  continuous  at  the  wave  front.  Then 
"i  -  “2  ’  0 
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By  differentiating  this  equation  along  the  wave  front  in  r -direction 
one  can  write 

(u{-  up  +  i  (ux-  u2)  =  0 

where  (  ) '  denotes  partial  differentiation  with  respect  to  r  and  (  ) 
denotes  partial  differentiation  with  respect  to  t  .  Then,  from  Eq.  (5.^1) 
one  has 


2 

c 


1 

P 


°r2~  °rl 
er2  “ €rl 


(5^2) 


if  both  the  sides  (l)  and  (2)  are  in  the  locked  region  and  the  locking 
is  nonideal.  One  can  derive  the  following  value  of  the  wave  speed  from 
Eq. (5.42) ,  the  stress-strain  relationship  (5*35)  and  the  displacement 
continuity  relationship,  namely: 

C  . 

</  3p 

This  is  the  same  value  as  that  obtained  in  Eq. (5-39).  Further,  as  in 
the  case  of  one -dimensional  wave  propagation,  when  region  (2)  is  elastic 
and  region  (l)  is  locked,  it  can  be  shown  that  the  loading  waves  catch 
up  with  each  other  and  move  as  a  shock  wave  with  finite  amplitude,  while 
the  unloading  waves  fall  apart.  Then  that  wave  which  brings  the  material 
to  a  state  of  incipient  unloading  propagates  with  the  speed  just  given, 
and  it  is  followed  by  an  elastic  wave  of  speed  ^e  +  2(^p'  . 


(ii)  Infinite  Medium  with  Pressure  Suddenly  Applied 
at  the  Edge  of  a  Spherical  Cavity 

In  this  section  we  shall  consider  an  infinite  medium  (made 
of  a  volumetrically  locking  material)  with  a  spherical  cavity  of  radius 
a  as  indicated  in  Fig.  5.10.  The  stresses,  strains  and  displace¬ 
ments  in  the  medium  due  to  a  pressure  p  suddenly  applied  at  the  cavity 
surface  will  be  studied.  Further,  it  is  assumed  that  the  pressure  p 
is  maintained  at  the  cavity  surface  for  all  times,  i.e.. 
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r=a 


p  for  t  >  0  . 


When  the  pressure  p  is  suddenly  applied  at  the  cavity  surface 
a  wave  starts  spreading  into  the  medium.  At  a  certain  time  t  this 
wave  front  will  he  at  a  certain  radius  r  =  k(t)  .  Then,  the  region 
a  §  r  £  k(t)  will  be  a  volumetrically  locked  region  if 

(a)  the  pressure  p  is  of  sufficient  magnitude  to  cause 
volumetric  locking, 

(b)  the  strains  in  the  region  0  S  r  S  k(t)  are  nowhere 
of  a  magnitude  which  would  cause  unlocking. 
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The  preceding  requirements  will  be  checked  later. 


The  stresses,  strains  and  displacements  in  the  region  0  i  r  S  k(t) 
are  governed  by  the  following  equations;  (equation  of  motion,  distor- 
tional  elasticity,  locking  condition  and  kinematics) 

h5r  ^  .  V  de  du 

(5.43) 

(5.44) 

?r  +  2?e  -  *  ei 

(5-45) 

€0 

The  locking  is  assumed  to  be  ideal.  Further,  it  is  assumed  to  take 
place  in  compression  only.  Then  the  locking  condition  (5.45)  holds 
throughout  the  locked  region.  Eqs.(5.45)  and  (5.46)  then  yield  the 
following  differential  equation  for  u . 


Su 

5? 


u 

r 


(5.46) 


dtt  2u  _ 

5 r+T  +  U  =  0 


with  the  general  solution 


1  = 


where  f(t)  is  an  arbitrary  function  of  time.  Then, 


2f(t)  _£ 

3  *3 

r 


(5-47) 


(5.48  a) 
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(5.48b) 


ee  “ 


f(t)  _  h 


Then,  from  distortional  elasticity  one  obtains  the  following  relationship 

0r  -  ae  =  -  6g  ^  (5.49) 

r 


The  equilibrium  equation  and  this  equation  then  yield  the  following 
differential  equation  for  a 


12G 


sMfom 


A  simple  integration  yields  the  following  expression  for  the  stress: 

a  =  -  4g  -  p  +  sfa)  (5.49a) 

r  r  r 

Eqs . (5.47),  (5-48)  and  (5.49)  constitute  the  general  solution  in  the 
region  0  £  r  £  k(t)  .  The  solution  is  determined  if  the  functions 
f (t) ,  g(t)  and  k(t)  are  determined.  Therefore,  the  investigation 
in  the  next  section  will  be  concerned  with  determining  these  functions. 

(iii)  Functions  f (t) ,  g(t)  and  k(t) 

For  t  >  0  the  boundary  condition  at  the  cavity  surface 

r  =  a  is 

or  -  -  p  (5.50) 

At  the  wave  front  the  displacements  are  continuous.  Further,  it 
is  assumed  that,  for  the  time  t  <  0  ,  the  infinite  medium  is  at  rest 
and  is  stress  free.  The  pressure  p  is  applied  at  time  t  =  0  .  Then 
the  region  ahead  of  the  wave  front  r  *  k(t)  is  at  rest  and  is  stress 
free,  i.e.,  at  r  *  k(t) 


u  -  0 


(5.51) 
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and  hence 


0 


Further,  the  region  0  <  r  <  k(t)  is  in  a  state  of  ideal  locking  and, 
therefore  at  r  =  k(t) 


whence 


Hien,  from  Eq.  (5*^2) 


(5-52) 


(5*53) 


Hie  conditions  (5*50),  ( 5 • 5l)  and  (5-53)  will  be  used  to  determine 
the  functions  f(t),  g(t),  and  k(t).  From  Eqs.  (5*50)  and  (5*^9),  one 
has 


a 

Then 

«r  -  -  ijj  *  pf  p 

From  Eqs.  (5-5l)  and  (5*^7)  one  obtains 


(5* 53a) 

(5-5^) 

(5.56a) 

(5.56b) 
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e  * 
r 


3  3  t 


ii ti 

3 


(5.56c) 


(5-56d) 


0 

r 


4G 

3 


( 5 • 56e ) 


Then  from  Eqs.  (5*53)  and  (5.56e)  the  following  differential  equation 
for  k(t)  can  be  obtained: 


k 


0 


(5-57) 


If  the  function  k(t)  is  obtained  by  integrating  this  equation,  the 
functions  f(t)  and  g(t)  can  be  determined  from  Eqs.  (5*56a)  and  (5>53a). 
The  initial  conditions  for  the  differential  equation  (5*57)  are 


k(0)  =  a 


-P€/^2)t=0  =  [or] 


t=0 , r=a 


=  -P 


(5-58) 


Before  proceeding  to  integrate  the  differential  equation,  the  condition 
under  which  the  region  0  <  r  <  k(t)  could  be  a  locked  region  will  be 
investigated  in  the  next  sub-section. 

(iv)  The  Conditions  Under  Which  the  Locked  Region  0  <  r  <  k(t)  Exists 

The  region  0  <  r  <  k(t)  can  be  a  locked  region  if  the  stresses 
satisfy  the  inequality 

-(cr  +  2oQ)  >  3K€/  (5.59) 

where  3K  is  the  bulk  modulus 


K  =  3(1  -  2v ) 
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From  the  distortional  elasticity 
Then 


Sr  -  M<«  -  €,)  ■  c  + 


Oq  is  related  to 

6Gf(t) 

3 


a 


r 


(see  Eq.  5-*+9). 


Ihen  the  inequality  (5*59)  can  be  written  in  the  following  way  (see 
Eq.  5-5*0 


3p 


3pi-  (A  .  i)  -  120  %  >  3Ke, 


(5-6o) 


In  order  that  the  locked  region  a  <  r  <  k(t)  does  exist,  the  inequality 
must  hold  at  least  for  r  =  k  =  a,  t  =  0,  i.e-, 


P  > 


3K  +  4G 


(5-60a) 


for  v  =  0.25,  3K  =  56-  Then  this  inequality  becomes  p/G©f  >  3.  Also, 
from  the  inequality  (5-60),  one-can  see  that  the  material  will  unlock 
if,  at  any  place, 

3p  -  3pf  (i  -  =  K€f  (5.60b) 

This  will  occur  first  where  r  is  maximum,  i.e.,  at  r  =  k(t). 

Further,  from  Eq.  ( 5 ■  *+2 )  and  the  incipient  locking  conditions,  the 
following  value  for  dk/dt  at  the  instant  of  unlocking  can  be  derived 


dk  I  °r2  "  °rl  1 

dt  =  If  €r2  *  6rl  ? 

Material  ahead  of  the  wave  is  at  rest  and  stressfree.  Then  a  „  =  e  0  —0. 

r2  r2 

Therefore 
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are  related  by  the  elastic 


At  the  instant  of  unlocking  and  e  ^ 

law.  Then 


d*  Ji  V6ri*2'ei)+S0en' 

df|?  «rl 

where  \  is  Lamp's  constant.  Also  at  r  =  k(t) 
e 


Then 


dk 

dt 


which  is  the  speed  of  the  elastic  irrotational  wave. 


(v)  Solution  of  the  Differential  Equation  ( 5 . 57 ) 

Equation  (5-57)  is  an  ordinary,  nonlinear  differential  equation 
for  k(t).  It  may  be  solved  by  numerical  integration.  However,  before 
any  numerical  procedure  can  be  started,  it  is  necessary  to  clarify  the 
behavior  of  the  solution  at  k  =  a.  This  will  be  done  first,  and  then 
Milne's  method  will  be  used  for  numerical  integration.  Also,  to  prepare 
the  equation  for  numerical  work  it  is  non-dimens ionalized  by  Introducing 
the  following  new  variables 


In  this  notation,  Eq.  (5*57)  reads  as  follows. 

q(q  -  i)  +  (2q  -  l)  +  (q3  -  l)  -  a2  =  o 

where 

=  r-2- 
7oti 


(5-61) 


The  initial  conditions  for  this  differential  equation  are 
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T  =  0 


(5-62) 


k.  . 

q  =  -=  1 


and  (see  Eqs.  5*53,  5?58,  and  5*62) 


T  =  0 


dq  5pa  1  dk 
dT  ol;adt=a 


(5.63) 


The  power  series  solution  is  assumed  in  the  following  form, 
q  =  1 


+  OUT  +  b.T2  +  b_T2  +  b.T*  +  bcT^ 


5 


(5-64) 


Because  this  power  series  will  be  used  to  evaluate  the  values  of  q  and 

its  derivatives  for  T  <  0.003,  it  is  reasonable  to  expect  that  the  first 

2  2 

few  terms  of  the  series  should  be  sufficient  to  evaluate  d  q/dt  to 
the  accuracy  of  three  decimal  places.  Therefore,  the  coefficients  of 
q6  and  higher  powers  are  neglected  in  the  series. 

Then  by  keeping  only  terms  up  to  the  fifth  and  lower  powers  of  T, 
one  can  evaluate  the  individual  terms  of  the  differential  equation  (5.61). 
By  equating  like  powers  of  T  the  coefficients  bg  ...  b^  can  be 
calculated. 

For  a  particular  value  of  p/Ge.^  =  4  the  coefficients  are  as  fol¬ 
lows  (From  Eq.  5 -60a ,  for  V  =  0.25,  one  can  see  that  the  value  of 
p/G€£  creates  a  locked  region.) 


»3-  •'» 

»U  »  -3-6 
b?  =  8.09 


The  expression  for  q  can  then  be  written  in  the  following  way 

q  =  1  +  ST-It^+VFt3-  S.ST*  +  8.09  T5 
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This  power  series  has  been  used  to  start  the  numerical  integration. 

The  typical  numerical  solutions  for  p/G =  4  is  shown  in  the  follow¬ 
ing  table.  Unlocking  takes  place  at  q.  =  1.142  and  T  =  0.088.  The 
stresses,  strains  and  displacements  have  been  plotted  in  Fig.  5-12. 

Fig.  5 • 11  shows  the  position  of  wave  front  and  wave  velocity. 


r 

4 

4 

9 

0.0 

-3.0 

1.73206 

1.0 

0.01 

-2.899 

1.70257 

1.01717 

0.02 

-2.809 

1.67404 

1.03405 

0.03 

-2.714 

1.64649 

1.05066 

0.04 

-2.650 

1.61962 

1.06700 

0.05 

-2.574 

1.59353 

1.08305 

0.06 

-2.505 

1.56813 

1.09887 

0.07 

-2.441 

1.54341 

1. 11441 

0.08 

-2.381 

1.51930 

1.12974 

unlocks* 

0.09 

-2.328 

1.49577 

1.14480 

* 

Unlocks  at 

q.  =  1.142 

4  =  1.5 

In  the  next  sub-section  the  problem  after  unlocking  will  be 

studied. 
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-  l4l  - 


m|£l 


-  1^3  - 


Fig. 


It  has  been  shown  (see  p.  140  )  that  unlocking  begins  at  the  wave 

front  at  point  A  (Fig.  5*13)  where  k(t^)  =  r^«  When  the  wave  has 

reached  this  point,  its  velocity  has  the  value  c  of  elastic  waves, 

a b  everywhere  on  the  wave  front  u  =  0,  and  the  stress  at  A  is 
2 

=  pc  From  here  on  the  wave  continues  as  an  elastic  wave  going 
outward,  as  shown  in  Fig.  5-12. 

The  region  behind  the  elastic  wave  may  contain  a  locked  region  (2) 
and  an  elastic  region  (3),  as  indicated  in  the  figure-  The  wave  front 
AB  is  a  straight  line,  its  slope  corresponds  to  the  elastic  wave  speed 
c.  The  problem  is  now  to  obtain  the  stresses,  strains  and  displacements 
in  the  elastic  and  locked  regions  such  that  they  match  on  the  interface 
AD  which  is  still  unknown.  Further,  the  solution  must  satisfy  the 
appropriate  boundary  and  initial  conditions. 

The  general  solution  in  the  locked  region  is  still  given  by 
Eqs.  (5.^7),  (5>^)  and  (5.49).  The  general  solution  in  the  elastic  :• 
region  will  be  derived  in  the  next  section. 

(vii)  General  Solution  in  the  Elastic  Region 

The  dynamic  response  in  the  elastic  region  is  governed  by 
the  following  equations: 
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(a)  Equation  of  Motion  (aee  Eq. (5.34) 
a  -  a  v2 

2^— &  » 


Sar 

ST"  + 


Ste 


(b)  Stress -Strain  Relationship 


\  (e  +  2ea)  +  2Ge 
e'  r  Q'  r 


=  Xe(€r+2e©)  +  2G€6 


where  X.  is  Lamp’s  constant  and  2G  Is  the  shear  modulus  of  the 
e 

material . 


(c)  Kinematics 


€ 

r 


Su  u 

Sr  ’  6  =  r 


From  these  equations  the  following  differential  equation  in  u  can  be 
derived 


S2u  2  Su 

72  +  75?  - 


p  d2u 

\  +  2G  St2 


As  explained  in  p.128,  this  equation  can  be  rewritten  in  the  following 
form  by  introducing  the  potential  function  i(r  such  that  u  »  di|r/dr  : 

7  a  [L  /2  =  a  aft 

5r  L2  5r  If  Sr  J  ^  3t2 

where 

/^T2G 

c  =  Y  p 

Integration  of  this  equation  yields 
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7^89 


where  is  a  constant.  This  equation  can  also  he  written  as 


(5.66) 


The  general  solution  of  this  equation  is 

C  r^ 

r\|r  =  f^r+ct)  +f2(r-ct)  + 


(5.67) 


Once  is  known,  the  displacement  u  and  hence  the  stresses  and 
strains  can  he  calculated. 


We  are  seeking  solution  in  a  region  which  extends  to  infinity. 
Thus,  there  are  no  incoming  waves.  Further  the  elastic  waves  start 
at  time  t  =  t^  ,  and  at  radius  r  =  r^  .  Therefore,  the  potential 
function  i(r  in  the  elastic  region  can  he  written  in  the  following 
form: 


?f3  * 


(5.68) 


In  the  further  analysis  the  ahhreviation  t  =  t  - 1^  will  he  used. 
Then 


(5.69) 


where  (  ) '  denotes  differentiation  with  respect  to  the  argument  of 
the  function. 


(viii)  Boundary  and  Interface  Conditions 

The  general  solutions  in  the  locked  and  elastic  regions  are 
given  hy  Eqs.(5.47)  to  (5.49)  and  (5.69).  They  contain  arbitrary  func¬ 
tions  f(t),  g(t)  and  fj  .  These  functions  and  the  unknown  radius 
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r  «  r2 (t)  of  the  Interface  can  be  determined  from  the  following 
conditions. 

(a)  At  the  wave  front  the  displacement  is  zero;  i.e. , 
at  r  =  r^  +  ct  : 

u  *  0 


(5.70) 


(b)  At  the  surface  of  the  cavity,  the  radial  stress  is  -  p  ; 
i.e. ,  at  r  =  a  : 


(5.71) 


(c)  At  the  interface  the  radial  stress,  the  particle  velocity 
and  the  displacement  are  continuous;  i.e.,  at  r  =  r^(t)  : 


u  =  u 


(5.72a) 


a  -  a 
r  r 


(5.72b) 


L 

u  *  u 


(5.72c) 


Eqs. (5.72b)  and  (5.72c)  are  not  independent  of  each  other  because  the 
continuity  of  the  radial  stress  implies  the  continuity  of  the  particle 
velocity  across  the  interface. 

(d)  Further,  the  stresses  in  the  locked  region  are  in  a  state  of 
incipient  locking  at  the  interface;  i.e.,  at  r  =  r2(t) 


°r  +  2a0  -  -  3Kei 


(5.73) 


From  this  equation  and  the  condition  that  the  strains  in  the  locked 
region  satisfy  the  locking  condition  (5.^5),  we  have  at  r  =  rg(t)  : 


of  +  2o0  =  3K(er+2e0) 


(5.7^) 


Also,  from  distortions!  elasticity  in  the  locked  region  at  r  =  r2(t) 


5r-5e-  2O(Jr-i0) 


(5.75) 
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From  the  fact  that  the  elastic  stresses  and  strains  satisfy  Hooke's  law 
we  have  at  r  =  rg(t) 

a  +  2  oa  =  3KU  +2€fl)  (5.76a) 

T  o  TO 

°T-°e  -  2C(«r-<e)  (5-76b) 

From  Eqs. (5.75)  and  (5.76b)  at  r  =  rg(t) 

(5r-or)  -  (5g-o0)  =  20  Qir-  €r)  -  (ee-«e)] 

From  Eqs. (5.72a)  and  (5. 72b)  this  equation  becomes  at  r  *  rg(t) 

Oq  -  Oq  *  -  2G(er-  €y)  (5*77) 

Similarly  from  Eqs. (5-7*0  and  (5 -76a)  at  r  =  rg(t) 


°e  ‘  °e  =  2  ^€r  "  €r^ 


(5.78) 


From  (5.77)  and  (5-78)  one  obtains  at  r  =  rg(t) 
0 


C  -  € 

r  r 


and  then  =  c_  because  of  continuity  of  the  radial  displacement  u 

O  U 

at  the  interface.  Then  at  r  =  rg(t) 

€r  +  2ee  •  *r  +  2ie  =  ’  ei  (5*79) 

This  equation  can  be  used  in  place  of  one  of  the  equations  (5-72a,b,c) 
or  (5-73). 

From  Eq.(5.69)  one  can  write  the  following  expressions  for  strains 
and  velocity  in  the  elastic  region: 
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(5.80a) 

(5.80b) 

(5.80c) 


Then,  from  Eqs. (5-79) ,  (5.80a,b) 


TO) 


■  •  ■  ciV2 


(5.81) 


This  is  a  useful  relation. 


Now  let  us  consider  the  remaining  boundary  conditions.  From  Eqs. 
(5.69)  and  (5.70) 

1  *s«>)  Cl 

-  — — =— -r  fi(o) - - - s  +  ^t^  +  ct)  -  0 

=  (r1  +  c«  5'  (,  te if  51 


1*6*  ^ 

f‘(0)  -  0 

f3(0)  -  0  (5.81a) 

C!  -  0 

Then  Eq. (5.8l)  becomes 
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(5.82) 


"  eicr2 


From  (5.72a),  (5.89)»  and  (5-47)  one  has 

-  -  £  'J  ({  -^) 


r“(t) 


(5.83) 


Similarly  from  (5. 72c),  (5.80c),  (5*47)  and  (5.8la) 

m  .  -if.  ({.ip) 

_2  r0c  3  \  c  /  2  3  \  c  / 


(5.84) 


From  (5.71)  and  (5.49a) 


g(t)  =  4g  +  p  -  p  (5.85) 

O  €1 

a 

Then,  as  before, 

5r  =  4Gf(t)(4-^5 )  +  p3?(t)  (j  -  j)  -  p  (5.86) 

This  equation,  distortional  elasticity  (5*44),  Eq.  (5.48)  and  the  boundary 
condition  (5-73)  result  in  the  following  equation  (more  simply  from 
(5.60b)): 

e-80 

Eqs. (5.83),  (5.84)  and  (5-87)  are  three  equations  for  three  unknowns 
f  ,f  and  r2  .  However,  the  use  of  the  equation  (5-8l)  simplifies 
the  analysis. 
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From  Eqs. (5.82)  and  (5.84)  one  obtains 


/  r2(t)  -  r.\  o 

f3  [*-— Z - -]  -  *£c4  ~  t)  (5.88) 

Now  from  this  equation  and  Eq.(5.83)  one  has  the  following  expression 
for  f3 


f 


3 


£  €  r3 
3  e£~2 


(5.89) 


Eq. (5-88)  is  a  differential  equation  in  t  which  is  the  independent 
variable.  If  both  sides  of  this  equation  are  differentiated  with 
respect  to  t  ,  one  obtains  the  following  equation 


C.ifa 

In  «  il 

L  cT0 

But,  from  Eq.(5.82) 

f5  V  «  I 

=  -  v 

dr2  ..  _ 


(5-90) 


Then,  from  Eq. (5*90 )  and  this  equation,  one  can  write  the  following 
expression  for  f(t)  : 


f(t)  =  e^c  rg 


(h ?  +  1)  <5-91) 

If  this  expression  for  f(t)  is  substituted  in  Eq. (5-87)  one  obtains 


.  f  (€) 

“TT  + 

a 


•vS&SMH) 


p  -  Ke, 


i.e. 


f(t)  ■  (p  -K€j)  -  pc^c‘ 


(5,92) 
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Then  the  following  expressions  for  i  and  f  can  he  obtained  by 
differentiating  this  equation  with  respect  to  t 


f(t) 


-  PV 


2  a 


(5.93a) 


f(t) 


pe/&3[  2  dr2  dS  1  dr2  r2  dS 

~TG  [_ac  dl2  +  ac  &i2  +  ac  d£3 


,,3  .,2  _ 

=  di3  8  at2  J 


(5.93b) 


By  equating  the  two  different  expressions,  (5.91)  and  (5. 93b),  for 
f(t)  the  following  differential  equation  for  r2(t)  can  be  obtained 


(5.9*0 

If  this  equation  is  solved  for  r2(t),  f(t)  can  be  obtained  from  the 
Eq.(5»92).  Then,  g(t)  and  f3  can  be  obtained  from  Eqs. (5.85)  and 
(5.89).  The  displacements,  stresses,  and  strains  can  then  be  determined 
in  elastic  and  locked  regions. 


(ix)  Solution  of  the  Differential  Equation  (5.9*0 

The  differential  equation  (5-9*+)  can  be  non -dimens ionalized 
by  introducing  the  following  variables 
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> 


l.e. 


T  -  0  ,  t  *  t^  or  T  ■ 


Then  the  equation  can  be  written  in  the  following  form  if  the  value  of 
the  Poisson's  ratio  is  assumed  to  be  equal  to  0.25,  i.e.,  3K  *  5G  and 
(see  Eq. (5.66a) ) 

/3£  /3K  +  4G  3p  3 

C  V^G  =  V  3p  k5  =  2 


d\ 


d  R  d2R  /  dR  .  dlU 


dR2 

3  ——  + 
dT 


**2 


+  2  R2 


*  0 


(5-95) 


The  initial  conditions  for  the  differential  equation  can  be  obtained 
from  the  conditions  at  the  time  of  initial  unlocking  at  the  radius 
r  =  ,  i.e.,  r  =  r^  and  t  =  t^  or  t  =  0  ,  in  the  x  -  t  diagram 

(Fig.  5- 13). 

For  example,  one  can  consider  the  case  when  p/Ge^  =  4  .  Then, 
from  Fig.  5  •  Ha 

R2(0)  »  1.-142 

Now,  the  values  of  dR(0)/dT  and  d2R(0)/dT2  can  be  obtained  from 
Eqs. (5.92)  and  (5. 93a)  in  the  non-dimensionalized  form,  i.e.. 


Again,  the  values  of  dRg/dT  and  d2Rg/ dT2  in  the  case  when  p/Ge^  =  4 
are 
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T  «  0 


-  0.2676 


25.229 


> 


Then,  Eq.(5*86)  can  be  integrated  numerically.  Milne's  method  is 

used  for  numerical  integration.  The  typical  numerical  solution  has 

been  illustrated  in  the  following  table  for  the  value  of  p/Ge,  =  4  : 

t 


- 

i \ 

d\ 

T 

R2 

dT 

df2 

1? 

0 

1.142 

-0.2676 

-25.229 

94.133 

0.005 

1.1404 

-0.3930 

-24.92 

29.98 

0.01 

1.1381 

-0.5173 

-24.93 

-33.66 

0.015 

1.1352 

-0.6428 

-25.26 

-101.63 

0.020 

1.1317 

-0.7705 

-25.96 

-178.76 

0.025 

1.1275 

-0.9031 

-27.07 

-272.58 

0.030 

1.1227 

-1.0421 

-28.73 

-393.36 

0.035 

1.1171 

-1.1915 

-31.08 

-559.25 

0.040 

1.1108 

-1.3546 

-34.45 

-801.84 

0.045 

1.1035 

-1.5384 

-39.31 

-1180.07 

*  unlocking  along  GH  (Fig.5.l4) 

The  variations  of  rg  ,  the  stresses,  strains  and  displacements  in 
the  region  are  illustrated  in  the  Figs. 5-14  and  5.15*  It  can  be  seen 
that  unlocking  takes  place  throughout  the  region  GH  (Fig. 5*14)  at  time 
T  *  0.044  ,  or  T  «  1.32  .  The  solution  for  the  problem  must  be  vorked 
separately  for  time  T  k  1.32  . 


-  155  - 


Fig.  5 -15a 
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